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METEOROID HAZARD EVALUATION FOR SIMPLE STRUCTURES 
WITH VARIOUS ORIENTATIONS 
byC. D. Miller 
Lewis Research Center 

SUMAAARY 

An analysis was performed on the effects of orientation of planes and cylinders upon 
collision hazard due to meteoroids of mass within the photographic range at a position 
in Earth's orbit, but far from Earth. Use was made of previously developed weighting 
factors to correct for biases in the photography of meteors. These biases were caused 
by meteoroid mass and velocity and the direction in space from which a meteoroid 
comes. The analysis was based upon data for photographic meteors supplied by the 
Smithsonian Astrophysical Observatory. 

A damage criterion supplied by Ames Research Center was used; it contained as 
factors the normal component of impact velocity to the 0. 875 power and particle mass to 
the 19/54 power. 

The analysis showed a substantial advantage in (1) orienting a plane sheet of armor 
normal to the apex of Earth movement, (2) orienting a cylindrical tube in a direction 
with one end of its axis pointing toward the Sun, and (3) arrangement of a panel of tubes, 
with one end of the axis of each pointing toward the Sun, all in a plane normal to the apex 
of Earth movement. However, the relative advantages of orientation would rapidly 
change with changes in the ratio of the exponents of mass and velocity in the damage 
criterion used, possibly even to the extent of reversing. 

Bimodal log-normal equations were obtained and integrated for normal components 
of impact velocity on planes, cylinders, and spheres with three principal orientations of 
the planes and cylinders. Necessary armor thicknesses for a desired probability of no 
destructive impact were determined with a double integration of particle mass and ve- 
locity so that all possible combinations of mass and velocity were taken into account. 

Implications are found and discussed to the effect that a tentative description of the 
total meteoroid population of the solar system can be derived, from various character- 
istics of the sample existing in Earth's orbit. By similar methods, the overall meteor- 
oid hazard could be estimated for an interplanetary mission. 
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INTRODUCTION 

Meteoroids, small particles moving in orbit about the Sun, constitute a hazard to 
various components of space vehicles. Space radiators that might be required for a long 
mission are particularly vulnerable to destructive impact. Such radiators, used for re- 
jection of waste heat from powerplants,, would need to be large because their effective- 
ness would be based entirely on radiation of heat to space, not including any effect of 
conduction and convection as with land-based powerplants. 

A necessary step in calculation of the hazard from meteoroids has been determin- 
ation of the flux rate relative to mass and of the velocity distribution for a given mass, 
because the greater the mass or the velocity of a meteoroid impacting upon a metal 
surface, the greater the damage. 

Meteoroids that have relatively high probability of destructive impact upon space 
radiators are those within the so-called photographic range of mass. When such a me- 
teoroid encounters the atmosphere of Earth, ablation or rapid erosion occurs with much 
luminous radiation. The result is a streak, known as a meteor, having sufficient lumi- 
nous energy to permit photography from the ground. 

As explained in reference 1, prior to the work reported therein, other investigators 
had advanced meteor theory to a point such that relative values of mass and absolute 
values of velocity, orbital parameters, and other characteristics of the particle that 
caused a given meteor could be well determined. (Absolute values of mass are still 
somewhat imcertain. ) Difficulties still existed because of important biasing effects in 
the photographic process caused by (1) meteoroid mass, (2) meteoroid velocity relative 
to Earth's atmosphere, and (3) zenith angle (angle of the meteor path through the at- 
mosphere to the zenith). An additional spatial bias existed because of the necessary lo- 
cation of the cameras in a fixed position (New Mexico) rather than in positions randomly 
varied over Earth's surface, the necessity of operation of the cameras exclusively at 
night, and the operation of the cameras principally during the later hours of the night. 

Wei^ting factors to correct the various biases were developed and well confirmed 
in the work reported in references 1 to 5. With use of these wei^ting factors, velocity 
distributions without regard to direction of impact were determined both relative to a 
hypothetical gravity-free Earth and relative to the real Earth. Log-normal equations 
were developed representing both velocity distributions closely. 

Also, the equation for meteoroid influx relative to mass from reference 6 was well 
confirmed throughout the photographic range. A double integration for meteoroid damage 
was performed with use of (1) any specified frequency fvmction for velocity (velocity dis- 
tribution), (2) the mass influx equation, and (3) a generalized equation for use with any 
damage criterion provided such criterion is proportional to the product of some power of 
mass and some power of impact velocity. 
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In the work reported here, the weighting factors were used in determining the dis- 
tribution of normal components of impact velocity on spheres and on planes and cylin- 
drical surfaces oriented in the three principal manners relative to the ecliptic plane and 
the apex of Earth movement (the momentary direction of Earth’s movement). Log- 
normal equations for velocity distribution and sample calculations of necessary armor 
thickness are included. As in the earlier work, use was made of data from reference 7, 
which were provided on punched cards by Smithsonian Astrophysical Observatory. 

Attention will be given to the bearing of the results of the analysis on the possibility 
that the meteoroid population at any desired position in the solar system may be tenta- 
tively deduced. 


METHOD OF PROCEDURE 

Distributions of normal component of impact velocity in the form of histograms were 
constructed for surfaces of structures traveling in Earth’s orbit about the Sun but 
assumed to be far enough from Earth that impacting meteoroids would be negligibly 
affected by Earth’s gravity. Three basic structural forms were included: (1) a sphere, 
(2) a plane, and (3) a cylinder having infinite ratio of length to diameter. For the plane 
and the cylinder three orientations were included as is illustrated in sketch (a) below. 
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An axis normal to the plane, or the axis of the cylinder, was alined with positive direc- 
tion as follows: (1) pointing in the direction of Earth movement in its orbit (i orienta- 
tion), (2) pointing in a direction approximately away from the Sun, but exactly perpen- 
dicular to the i orientation and within the plane of the ecliptic (j orientation), and 
(3) pointing in a direction normal to the ecliptic plane southward (k orientation). A de- 
scription of the method of computing data on distributions of normal velocity components 
for use in constructing the histograms just discussed, with reference to several appen- 
dixes, will constitute the first major section of this presentation. 

Log-normal equations were developed to represent at least approximately the dis- 
tribution of normal velocity component shown by each histogram. A description and 
discussion of these equations will be a second principal section of the presentation. 

A third section will deal with computation of a concentration factor <p for each 

o 

orientation of a plane or a cylinder, representing the frequency of impacts as compared 
with the expected frequency for a similar structure randomly oriented. 

The fourth principal section will present sample computations of necessary armor 
thickness for the structures and orientations treated. Quantitatively, these computa- 
tions will not be intended as a definitive result. They will, however, demonstrate a 
method that is generally applicable with a wide variety of damage criteria that might be 
selected and they will provide a comparison of the effectiveness of the various orienta- 
tions of the structures studied. The computations will use the log-normal distribution 
equations and the values of cp^ in the manner described in reference 4. The sample 

calculations required integration of the log-normal velocity distributions to allow con- 

th 

struction of curves representing the average value of the n^" power of the normal com- 
ponent of impact velocity (where is the normal component of impact ve- 

locity and n is any arbitrary positive number). 

The fifth principal section, and the last preceding a general discussion of results, 
will be concerned with the possible bearing of the study on the general meteoroid pop- 
ulation of the solar system. 


COMPUTATIONS FOR USE IN CONSTRUCTION OF HISTOGRAMS 
OF NORAAAL VELOCITY COMPONENTS 

Distributions of normal components of impact velocity were obtained as weighted 
counts of sporadic meteors within each of 80 velocity classes with use of a computer 
program. Each class had a central value (class mark) of an integer plus one-half, from 
1/2 to 79^ kilometers per second (km/sec). Thus, each velocity class extended through- 
out a range of 1 km/sec. The sporadic meteors used for the weighted count were a por- 
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tion of those for which data were published in reference 7. They were chosen in a man- 
ner that will be explained later. For explanation of the method of making the weighted 
count, it is convenient to use the concept of an adjusted unit count, which will now be 
explained. 

Adjusted Unit Count 

For determination of an adjusted imit count, the weighting factor developed in ref- 
erences 1 to 5 was used. 


(P^ = (cos F(Z 


R'av <= 


3.87 




( 1 ) 


As was explained in those references, the correct distributions of velocity or various 
other parameters should be obtained if each meteor is counted, not as one meteor, but 
as meteors. That is, cp^ may reasonably be termed an "adjusted unit count" for 
statistical purposes. 

In the right-hand side of equation (1), cos Zj^ is cosine of the zenith angle, or angle 
between the path of the meteor and the zenith (CZj^ of ref. 7). The symbol v^, as in 
reference 7, represents the velocity of the meteoric particle relative to Earth's atmos- 
phere before any deceleration by resistance of the atmosphere. The symbol F(Zd)„ 

xx 3^V 

represents a statistically expected value of a function of zenith angle, position of the 
meteor within the camera field of view, and azimuth of the meteor path, which was cal- 
culated and recorded for each meteor as part of the work described in reference 1. The 
value of F(Zj^)g^y is concerned with the density of the meteor trace on the photographic 
plates. The values recorded in the work of reference 1 were used here in equation (1). 
The value of cp is a spatial weighting factor, inverse values of which were calculated 

o 

in the work described in reference 5. It is concerned with the probability that a meteor 

would reach the atmosphere over the camera sites in New Mexico during the times at 

which the cameras were operated, in consideration of the direction in space from which 

it came and in consideration of its geocentric velocity. The probability <p for a meteor 

is not absolute, but is correctly related to the value of (p calculated for any other 

s 

meteor. (In general, the significance of any symbol used in this report is explained in 
appendix A. ) 

To obtain data for the desired histograms, the adjusted unit count was applied to the 
meteors that were selected for use. Orbital data given in reference 7 were used to de- 
rive components of meteor velocities normal to the surfaces of the structures that were 
under study. The adjusted imit counts, or fractions of them, corresponding to a given 
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class of normal velocity component were added to produce an adjusted total for that nor- 
mal velocity class. The types of meteors included in the analysis, the types of orbits 
treated, simplifying approximations made, the orbital parameters used, and certain 
needed basic constants will be discussed before the method of resolution of velocities is 
begun. 


Meteors Used for Determination of Velocity Distributions 

Of the sporadic meteors reported in reference 7, 1282 were used, as in reference 4. 
All rejections of sporadic meteors were for two reasons: 

(1) To avoid unrealistic values of the adjusted unit count cp^ of equation (1) meteors 
having a value of cos Zj^ less than 0. 2 were rejected. The manner of calculation of the 
inverse value of cp in the work of reference 5 included a provision to eliminate any 
statistical bias that might be introduced by these rejections. 

(2) Meteors were rejected if, at the time of impact, the meteoric particles were 
approaching the ecliptic plane from the south or if the components of their velocities 
approximately in the direction away from the sun and exactly perpendicular to the apex 
of Earth movement were positive. This provision, as explained in reference 4, confined 
the sample of meteoroids to a quadrant of the space surrounding Earth that could be 
completely covered; that is, at some periods during the operation of the cameras, 
meteors of any velocity could be photographed from any direction whatever lying within 
that quadrant. Restriction of meteors to that quadrant, as explained in reference 4, was 
based on an assumption of two planes of symmetry for approach of meteoroids toward 
Earth. 


Orbital Types and Simplifying Approximations 

A decision was made to use the equations for each of the three conic sections, 
ellipse, parabola, and hyperbola. Use of the three sections involved no decision on the 
controversial question as to whether parabolic or hyperbolic orbits actually exist. The 
data as presented in reference 7 indicate such orbits, and, even though such indication 
may well be due to inaccuracy of the data reduction, no rule could be set by which the 
given data could be warped to form elliptical orbits without the danger of introducing 
even greater inaccuracy. 

Hyperbolic orbits in reference 7 were indicated by negative values of a (semi-major 
axis, in case of ellipse) and the absence of values of q' (aphelion). The negative values 
were compatible with the other parameters given when made positive (for later purposes) 
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and used in the basic equation of a hyperbola 



( 2 ) 


Parabolic orbits were indicated by absence of values of both a and q’, with eccen- 
tricity e equal to 1. 00. The indication of a parabolic orbit is equivalent to an indica- 
tion that the heliocentric velocity Vjj is equal to the velocity of escape from the Sun. 

The values of Vjj given in reference 7 for the few cases of parabolic orbits were not 
identical, possibly because of difference of the Earth's distance from the Sun. In this 
analysis, a satisfactory degree of accuracy was possible with substantial simplification 
of the computations by treating the Earth's orbit as exactly circular with a radius of 
1 astronomical unit. Accordingly, for the parabolic orbits a calculated velocity of escape 
was substituted for the given value of Vjj in each case. For elliptical and hyperbolic 
orbits the value of Vjj (from ref. 7) was used. 


Orbital Parameters Used for Each Meteor 

Because of redundancy, parts only of the data presented for each meteor in refer- 
ence 7 were used. Certain checks were made on the consistency of the unused data with 
the data that were used. Details are presented in appendix B. 


Basic Constants 

Certain basic constants needed in the resolution of velocities are explained in 
appendix C. Derivations are included for the less obvious constants. 


Determination of Three Principal Components of Meteoroid 
Velocity Relative to Gravity-Free Earth 

As a first step in determining normal components of impact velocity of a meteor 
upon the structures under consideration, a vector representing heliocentric velocity 
was resolved into components extending in three directions: (1) the i direction, toward 
the apex of Earth movement in its orbit (positive for meteoroids in direct motion, nega- 
tive for meteoroids in retrograde motion), (2) the j direction, within the plane of the 
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ecliptic, approximately in the direction away from the Sun, but exactly perpendicular to 
the apex of Earth movement (negative for all meteoroids treated), and (3) the k direc- 
tion, normal to the ecliptic plane and southward (positive for all meteoroids treated). 
Such preliminary resolution of the velocity ^ is explained in detail in appendix D. 
Next, as also explained in appendix D, three similarly directed components of velocity 
were determined, but for the vector velocity v^ relative to the moving Earth rather 
than for the velocity ^ relative to the Sun. 


Determination of Velocity Components Normal to Surfaces of Structures from 
Principal Components of Velocity Relative to Earth 


Two of the components of namely, vT^ and from equations (D43) and 
(D44), were usable directly as normal components of impact velocities upon planes with 
the j and k orientations. Because of an assumption of symmetry relative to both the 
k and j planes, it was convenient and permissible to make the values of these com- 
ponents arbitrarily positive. The third component, vj^ from equation (D42), was usable 
directly as a normal component of impact velocity on a plane with i orientation. But 
here symmetry did not exist. If • i, with vj^ as given by equation (D42), was 
negative, impact occurred on the leading side of the plane. If that value was positive, 
impact was on the trailing side. 

Normal components of impact velocity on different parts of the surfaces of cylinders 
with i, j, and k orientations were readily determined with use of the vector components 

three vector components could have been used for 
determination of the normal component of impact velocity on various parts of a sphere. 
But, as a sphere presents the same shape, size, and relative orientation of vulnerable 
surface toward a meteoroid that approaches from any direction, the scalar value Vq as 
given in reference 7 was usable directly. 

The details of the determination of normal components of impact velocity on various 
parts of the surfaces of cylinders and spheres are explained in appendix E. 


Normal Velocity Distributions as Normalized Sums of Adjusted Unit Counts 

The circular symbols in figures 1 to 3 represent histograms. To avoid unnecessary 
confusion, the usual vertical bars of histograms are omitted. The circular symbols 
correspond to the tops of those bars. The curves, representing the general trend of the 
circular symbols, will be discussed in a later section. Effectively, the datum used for 
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plotting each circular point was a normalized sum of adjusted unit counts for all mete- 
oroids that would impact upon the pertinent structure (sphere, plane, or cylinder), with 
the pertinent orientation of that structure, with a normal component of impact velocity 
within the pertinent velocity class. The summation for each velocity class included only 
a part of each adjusted unit coimt . Such part of the adjusted imit count was equal to the 
total adjusted imit count (eq. (1)) multiplied by a relative probability that the particular 
meteoroid under consideration would impact the structure with normal velocity compo- 
nent within the pertinent velocity class. That relative probability consisted of two fac- 
tors. For definition of factor (1), consider a great number of hypothetical meteors, all 
exactly like the particular meteor imder consideration and all impacting the given struc- 
ture with the given orientation. Then factor (1) is the fraction of those hypothetical 
meteors that would impact with a normal component of impact velocity within the velocity 
class under consideration. Factor (2) is the vulnerable area presented to the given 
meteoroid by the given structure with the given orientation, measured as a projected 
area upon a plane normal to the path of the meteoroid. For planes, factor (1) for a 
given meteoroid was unity for one velocity class and zero for all others. For spheres 
and cylinders, factor (1) was nonzero for all velocity classes up to a maximum value of 
velocity. Sums of adjusted unit counts for the different velocity classes were finally 
normalized by dividing each of the sums by the total of all of them. Hence, factor (2) 
could be applied without reference to any particular standard size of structure so long as 
the dimensions of the structure were kept constant. The methods involved in obtaining 
the normal velocity distributions as normalized sums of adjusted unit counts are de- 
scribed in detail in appendix E. 


LOG-NORMAL EQUATIONS FOR DISTRIBUTION OF NORMAL COMPONENT OF 
IMPACT VELOCITY UPON PLANES, CYLINDERS, AND SPHERES 

In figures 1 to 3, curves have been constructed to represent unimodal or bimodal 
offset log-normal equations. Each triangular plotted symbol represents the result of a 
calculation, with use of the pertinent equation, for the percentage of all impacts having 
a normal component of impact velocity at the given value of the abscissa ±0. 5 km/sec. 
The curved line was drawn through the triangular points as guides. 

In each part of each figure, the equation represented is of the form 
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C.exp| 


+ 5 -) - [i. 
' norm v 


( 3 ) 


in which f(v ) is the fractional frequency of occurrence of a given value of v 

norm norm 

±0. 5 km/sec (multiplied by 100 for plotting in the figures). The offset 6. is arbitrarily 
equal to zero for i = 2. The constants Cj, and a^, as well as 6. with i = 1, may 
be varied independently to improve the fit of the resulting curves to the velocity histo- 
grams (circular symbols in figures). 

An equation of such form was fovind in reference 4 to provide a good representation 
of the gross impact velocities for both near-Earth and far -from -Earth (real Earth and 
gravity-free Earth) conditions. For figures 1 to 3, the constants of equation (3) were 
varied in a trial -and -error procedure to fit the equation as well as possible to the histo- 
grams. The resulting values of the constants are shown in table I. 

In four cases the fit of the curves to the histograms is so good, with due regard for 
the random vertical scatter of the data, that search for a better equation would not appear 
to be justified. Those cases are figure 1 (a, b, and c) and figure 2(c), with a minor ex- 
ception for figure 1(b) that will be discussed later. The four cases cover the leading and 
trailing sides of the i plane, either side of the j plane, and the k cylinder. In all 
other cases, the fit is obviously not perfect, though in each case it was the best that 
could be obtained with use of equation (3). Improvement could be effected with much 
labor by adding one or more terms to equation (3), that is, by varying i from 1 to 3, 


TABLE I. - VALUES OF CONSTANTS IN OFFSET BIMODAL LOG-NORMAL EQUATION TO PROVIDE BEST 
APPROXIMATION OF HISTOGRAM REPRESENTING NORMAL COMPONENT OF IMPACT 
VELOCITY ON STRUCTURES WITH VARIOUS ORIENTATIONS. 


Structure 

Coefficient in eq. (3) 

Offset 

Logarithmic modal value 

Logarithmic standard deviation 


Cl 

^2 

®i 

^^1 

1^2 

<^1 

<^2 

i plane (leading side) 

0. 1433 

1. 385xl0‘^ 

0 

1. 100 

3.90 

0.700 

0.200 

i plane (trailing side) 

. 1880 

0 

25.0 

3.420 

— 

.0692 


j plane 

.0883 

1, 365X10"^ 

5.0 

2. 560 

3.05 

.280 

.225 

k plane 

.0947 

2. 808X10"® 

2.0 

2.050 

3.29 

.470 

. 170 

i cylinder 

.0894 

8. 86X10"® 

: 2.0 

2.080 

3.04 

.450 

.225 

j cylinder 

. 1220 

2.910x10"“^ 

0 

1.600 

3.66 

. 560 

.270 

k cylinder 

.0899 

7. 53x10"® 

2.0 

2.050 

4.00 

. 503 

. 120 

Sphere 

.0919 

8. 13X10"® 

2.0 

2.010 

3.96 

. 510 

. 115 
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1 to 4, and so on. In all cases, however, the author believes the fit of the curves to the 
histograms, with only the two terms, is good enough to provide reasonable accuracy by 
use of the equations. 

In figure 1(b) four data points at 18^, 26|-, 28 ^, and 44^ km/sec, were ignored in 
the fitting of the curve because (1) these points lie on such a low level on the ordinate 
scale it was believed their inclusion would have negligible effect on the end result of the 
analysis, (2) they are so widely scattered that they do not constitute an adequate sample 
upon which to base a second mode (i = 2) for equation (3), and (3) the probability is high 
that the data points appear in the figure only because of error. 

The 18^-km/sec point was caused by meteor 3005 (p. 44, ref. 7); the 26|--km/sec 
point by meteor 11178 (p. 48); the 28^-km/sec point by meteors 4424 and 6977 (p. 67, 
and p. 42); and the 442-km/sec point by meteor 3204 (p. 47). Heliocentric velocities 
Vjj given for these five meteors are 66. 1, 83.0, 81.7, 61.2, and 77. 9 km/sec. These 
velocities correspond to hyperbolic orbits as is clear from comparison with the 42. 14- 
km/sec value (eq. (C6)) for a parabolic orbit, which is the boundary condition between 
elliptical orbits (particles permanently bound by the Sun) and hyperbolic orbits (particles 
arriving from or at least permanently departing for the space outside the solar system). 

If the heliocentric velocities reported for these meteors were erroneous by at least 
the amounts they exceed 42. 14 km /sec, the errors were great even for the graphical 
method used for the results reported in reference 7. Nevertheless R. E. McCrosky 
(of ref. 7) emphasized in a personal communication that a hyperbolic orbit for a meteor 
has never been confirmed when the data were subjected to precise reduction, and he 
thinks it unlikely that any of these five apparent hyperbolic orbits would be confirmed. 

Particles in hyperbolic orbits (correctly or not) affect all the histograms in figures 1 
to 3 and affect that in figure 1(b) in the regions below 18^ km/sec as well as at that 
velocity and above. But, except in figure 1(b), and only for the higher velocities there, 
no serious error results. An occasional multiplication of a correct nonzero frequency 
by a factor of 2 or 3 on the logarithmic scale of ordinates, because of a falsely high 
value of Vjj, is no more serious than the random scatter that is apparent in all parts of 
the histograms. 

CONCENTRATION FACTORS FOR FREQUENCIES OF METEOROID IMPACT 
UPON VARIOUSLY ORIENTED STRUCTURES 

In reference 4 a method was set forth for computing necessary armor thickness for 
a space mission to provide a specified low probability of avoidance of disastrous impact. 
Application of that method here will require the use of equation (3) with the values of 
constants shown in table I. It will also require, for each orientation of each structure, 
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TABLE n. - TOTAL IMPACTS ON PLANES, CYLINDERS, 
AND SPHERES AS PERCENTAGES OF EXPECTED 


IMPACTS ON RANDOMLY ORIENTED SURFACES. 


Structure and orientation 

Total impacts, 


<pg X 100 

i oriented plane (leading side) 

31. 10 

i oriented plane (trailing side) 

122. 30 

j oriented plane, either side alone 

121.07 

k oriented plane, either side alone 

91.26 

i oriented cylinder 

no. 30 

j oriented cylinder 

85.01 

k oriented cylinder 

101.78 

Sphere (in effect, randomly oriented surface) 

100.00 


a concentration factor <p as used earlier in reference 4. That factor may be defined 
as the ratio of the frequency of impacts upon the structure to the expected frequency of 
impacts upon one side of the same surface area when randomly oriented. 

Calculations of the concentration factor (p^ are explained in appendix F. The re- 
sults are tabulated in table n as percentages. 


SAMPLE COMPUTATIONS OF ARMOR THICKNESS FOR 
PROTECTION FROM METEOROIDS 

In reference 4 a method of computation that will be used here was presented and ex- 
plained for determining the necessary armor thickness to provide a desired high prob- 
ability of encountering no disastrous impact by a meteoroid during a space mission. 
Neither the sample computations in reference 4 nor those here are intended to provide a 
definitive result. In particular, the computations do not include any consideration of 
shower associated meteoroids, which would require further study. The computations 
are intended as a demonstration of a method that may be applied to various possible 
damage criteria as well as to a specific damage criterion used as an example. Also, the 
results here will serve as an illustration of the effect of structure orientation as govern- 
ed by a specific criterion and will indicate qualitatively how the effect of orientation 
might change with use of other damage criteria. 

As in reference 4, a space mission will be assumed with a duration of 1000 days, 
and the specific damage criterion used will be 
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= 0. 8162 
cr 


^-1/18 8/54 -1/2 19/54^,0. 875 


( 4 ) 


in which is minimum armor thickness (cm) at which penetration would occur, e^. is 
percent elongation of armor material, P„ is density (g/cm^) of particle, Po. in density 
(g/cm ) of armor material, m is mass (g) of particle, and is the normal com- 

ponent of impact velocity (km/sec). Equation (4) represents a formula provided by 
Richard H. Fish and James L. Summers of Ames Research Center in a personal com- 
munication. As in reference 4, stainless steel armor will be assumed with values of 
p^ equal to 8 grams per cubic centimeter and equal to 10 percent. Meteoroid den- 
sity pp will be taken as 0. 2 gram per cubic centimeter. 

In reference 4, the critical armor thickness necessary under the conditions spec- 
ified and with use of equation (4) was, in effect, found to be 


t^^ = 0. 200. 
cr 


3. 287X10"^^ TA 


(1. 34 e/X) 


X/1. 34 


/^g^norm 


(5) 


where T is the duration of the mission in seconds, A is vulnerable area in square 
meters, p is average number of destructive impacts expected per mission (very closely 
equal to the complement of the desired probability of no destructive impact when that 
desired probability is high), e is the exponent of ini ^ equation such as (4), and 

X is the exponent of m in the same equation. The symbol cp^ was used in reference 4 
with the same significance as here. The derivation of equation (5) involved a double 
integration relative to a well confirmed mass distribution and the velocity distribution 
shown by equation (3), so that impacts by meteors of all masses and all values of 
are involved in the solution. However, it was found that the result (eq. (5)) is usable 
for any velocity distribution. 

Equation (5) is simply a rearrangement of an equation that was derived in refer- 
ence 4, which was linear in a variable Z The critical thickness t ^ in equation (4) 
is a specific case for the general variable of reference 4. Now, let it be sup- 

posed that we might wish to use two or more expressions for t^^ similar to equation (4), 
each applicable only within a specific range of values of and each having its own 

values of e and X, that is, and Xj for the first expression, €2 and Xg for the 
second expression, and so on. Then, by the same method used in reference 4, we could 
easily derive a nonlinear equation in t^^ analogous to equation (5). With n equations 
like equation (4), for n different ranges of value, the nonlinear equation in t^^ 

would contain n terms in t^^, in which the symbol t^^ would be raised to n different 
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powers (integral or nonintegral) . The nonlinear equation would also include averages of 

n different powers of each average would be taken only over the pertinent 

range of v value. Such nonlinear equation could be solved for the critical value 

® norm 

ter approximate methods. Althou^ multiple roots would exist, physical reasons 
would also exist for excluding all roots but one. 

Substitution of numerical values in equation (5) yields 


( ' < ooo \0.2626 


( 6 ) 


In reference 4, methods were derived for integration of equation (3) under certain 
conditions to yield an average value v”„ . Under other conditions, recourse must be 

made to direct numerical integration. Use of those methods, including direct n umerical 
integration, with the constants shown in table I, produced the values of shown in 

figures 4 to 6, for planes, cylinders, and spheres, with the various orientations that 
have been discussed for planes and cylinders. 

Thicknesses of armor necessary for certain probabilities of no destructive impact 
and certain vulnerable areas were computed with use of equation (6) for planes and cyl- 

3 332 

inders, with the three principal orientations, and for spheres. Values of 

were taken from the integration results plotted in figures 4 to 6. The resulting armor 

thicknesses, with pertinent data, appear in table HI. 


TABLE m. - COMPUTED THICKNESSES OF ARMOR NECESSARY FOR PROTECTION FROM METEOROID 


IMPACT UNDER SEVERAL SETS OF CONDITIONS OVER A PERIOD OF 1000 DAYS 


Structure 

Orientation 

Vulnerable 

Approximate 

Average number 

Value of (pg 

Value of 

Critical thickness 



area 

probability 

of destructive 

from table V 

„3. 332 

(P^V 

by equation (6) 



.,2 

2 

of no 

impacts per 


g norm 





ft 

m 

destructive 

mission 



in. 

cm 





impact , 
percent 

(p in eq. (6)) 





Plane 

i (leading 

1000 

92.9 

99.5 

0.005 

0. 311 

7745 

0.209 

0. 530 



500 

46.5 

99.75 

.0025 

, 311 

7745 

.209 

.530 


i (trailing) 

1000 

92.9 

99.5 

.005 

1.223 

661.9 

. 109 

.278 



500 

46. 5 

99.75 

.0025 

1.223 

661.9 

. 109 

.278 


j 

2000 

185.8 

99.0 

.01 

1.211 

11 695 

.232 

.590 



1000 

92.9 

99.5 

.005 

1.211 

11 695 

.232 

. 590 


k 

2000 

185.8 

99.0 

.01 

.913 

6204 

. 197 

.500 



1000 

92.9 

99.5 

.005 

.913 

6204 

. 197 

. 500 

Cylinder 

i 

2000 

185.8 

99.0 

.01 

1. 103 

8752 

.215 

.547 


i 

2000 

185.8 

99.0 

.01 

.850 

5379 

. 190 

.481 


k 

2000 

185.8 

99.0 

.01 

1.018 

7243 

.205 

.521 

Sphere 


2000 

185.8 

99.0 

.01 

1.000 

6310 

. 198 

.502 
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IMPLICATIONS REGARDING TOTAL METEOROID POPULATION OF SOLAR SYSTEM 


A possibility exists that the meteoroid population of the solar system may possess 
interrelations necessitated by the effects of collisions between different meteoroids. 

These interrelations might be analogous to those existing between the molecules of a gas 
in a container. Although the collisions of meteoroids would of course be far from elastic, 
some changes in momentum components would occur upon collision and some kind of 
randomization of the motion should result. 

The frequency of collision for a given size of particle may be estimated with equa- 
tions presented in reference 4, of which equations (4) and (5) are specific forms. The 
general equations are 


Z = Cgm^^ 


( 7 ) 


and 


3.287x10'^^ TA 


^cr - ^31 




IX/1. 34 


( 8 ) 


where Z is any damage criterion, Z^^ is a critical value of Z, Cg is a constant, and 
V is an impact velocity of the applicable type. 

-3 

Now consider a particular spherical meteoroid of mass 10 gram and a density of 
0. 2 gram per cubic centimeter (see, e. g. , ref. 8 as source of density chosen). Assume 
the meteroid is moving about the Sun in the same orbit as the Earth. Let us estimate 
the time required for a 50-percent probability of at least one impact with another mete- 
oroid of the same or greater mass. 

The radius of the meteoroid would be 0. 106 centimeter. However, for impact with 

a similar meteoroid, its effective vulnerable area would be multiplied by four. In effect, 

its radius (i. e. , the total distance between centers of two meteoroids in contact) would 

be doubled. The vulnerable area, then would be 0. 1414 square centimeter, facing an 

approaching meteoroid from any direction. The corresponding total external vulnerable 

-5 

area of the sphere (A of eq. (8)) would then be 0. 5656 square centimeter (or 5. 656x10 

.3 

Now as any collision with a meteoroid of mass greater than 10 gram is to be con- 
sidered, the criterion Z in equation (7) is just the mass m. Hence, 
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C3=l 




X=1 > 


( 9 ) 


e = 0 J 


-3 

and in equation (8), = 10" ; that is, the critical mass m to be considered. For a 

50-percent probability of no impact within the specified mass range, according to the 
Poisson distribution, the value of p in equation (8) should be p = -In 0. 5 = 0. 6931. 

As the total vulnerable area is the surface of a sphere, in effect a randomly oriented 
surface, <p = 1. So, equation (8) becomes 

o 


10 


-3 


3. 287X10"^^ X 5. 656x10"^ T 


1/1.34 


0.6931 


( 10 ) 


and 


T = 3. 56X10^'^ sec = 1. 129xl0'^ yr (11) 

The time is indeed long, but very short in comparison with the age of the solar sys- 
tem. Moreover, collisions with smaller meteoroids would be many times more fre- 
quent. Therefore, it seems reasonable to assume that nearly all meteoroids have under- 
gone numerous collisions. Great heat would be generated, the resulting thermal explo- 
sions should release much gas, which presumably would become part of the solar wind 
and escape from the solar system. Any solid residue should be expected to be solidified 
froth (because of the expansion of hot bubbles of gas within it while in the molten state). 
Such condition would explain the low density of 0. 2 gram per cubic centimeter. But, 
most importantly, the residual particles should move in random directions relative to 
the average direction of motion of particles in Earth’s orbit. 

Figures 4 to 6 are interesting to examine relative to this collision hypothesis. It is 
obvious at once from figure 4 that the total flux of kinetic energy is at least approximately 
the same in both directions through the i plane, because the values of are 

the same for the two sides of this plane. (It was shown in refs. 2 and 3 that the velocity 
of meteoroids is statistically independent of mass. ) It is equally obvious that the flux of 
momentum (proportional to (^gV^^^^) is substantially greater for the trailing side of the 
i plane than for the leading side. 

Now consider the levels at n = 4 for curves that might be drawn through the four 
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sets of plotted values of ^g'^Qrm “ figure 4. For the leading side of the i plane, the 

level at n = 4 would be about 1. ixlO^ with = 0. 311. For the trailing side of the 

same p lane, the level would be about 2. 7x10^ with <p^ = 1. 223. The average value, 

for both sides of the i plane, then, would be the sum of these values of 

divided by the sum of the two values of cp^ or about 7. 3x10^. The le vel for the j plane 

from the figure is about 9. 4x10^ with cp„ = 1. 2107, so that the value v^ should be 

^4 ^g ^ norm ^ 

about 7. 8x10 . For the k plane, the value of <Pg^^Qj.j^ would be about 5. 6x10 with 

cp^ = 0. 9126, so that v^^^^ should be about 6. Ixl04. Now it is seen that the three 

values, 7. 3x10^, 7. 8xl04, and 6. 1x10^, are the same within a logarithmic range of 

about 0. 11. 

The three curves in figure 5 would yield about the same results. Althou^ these 
curves are for normal impact velocity on the surfaces of cylinders, each of them is 
equivalent to the average n " power of component of impact velocity on the family of 
planes passing through the axis of the pertinent cylinder. In figure 6, the curve for 
normal components of impact velocity on a sphere (equivalent to a randomly oriented 
plane) passes through approximately the same level at n = 4 as do each of the three 
curves of figure 5. 

Now, if we disregard particle mass because it is independent of velocity, the value 
v^orm proportional to the variance (square of the standard deviation) of kinetic 

energy. This variance relates to a zero value of kinetic energy (v = 0). In statistics, in 
a manner analogous to mechanics, the variance of v is known as the second moment of 
V. As in mechanics, the first moment becomes zero if recomputed in relation to its own 
value as a base. As with the second moment in mechanics (moment of inertia), the sec- 
ond moment can be corrected relative to a base equal to the value of the first moment 
simply by subtracting the square of the value of the first moment. 

Such a correctio n on the values of obtainable from figures 4 and 5, with use 

of the base values of v^^^^ that can be obtained in those figures, would produce very 
little change. Hence, the same level of for all the curves at n = 4 in figures 4 

to 6 is a direct demonstration that the standard deviation of kinetic energy is the same 
for the velocity components in all directions. This observation, together with the earlier 
mentioned observation from figure 4 (the equivalence of flux of kinetic energy from the 
two directions through the i plane), suggests that possibly the flux of kinetic energy 
should really be the same through all planes. The failure of figure 4 to show this fact at 
n = 2, then, could be due to the use of the wrong base value for v for the i and k 

planes. That is, possibly the j and k planes should be treated as moving relative to 
Earth. If they should be treated as moving southward and toward the Sun, then treating 
them as stationary relative to the Earth as has been done in this analysis would cause 
the flux of kinetic energy to appear too high through the j and k planes. This effect 
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would exist because this analysis was limited to meteors that approached the j plane 
from the antisolar direction and approached the k plane from the north. 

It is beyond the scope of this report to develop a kinetic theory for the highly in- 
elastic collision of hi^-velocity particles in a vacuum to show that the flux of kinetic 
energy should be invariant with orientation of the axes, even though the flux of momentum 
and flux of particles are both clearly not so. The purpose here is merely to examine 
the implication of that hypothesis, which has been suggested by the observations that 
have just been discussed. 

It is reasonable to assume that the average flux of kinetic energy of meteors through 
the Earth’s orbit has components southward through the ecliptic plane and inward toward 
the Sun. Such a flux would have to be compensated by outward and northward flux at 
other distances from the Sun. The suggestion of such a toroidal flux is not new (see, for 
example, the discussion of toroidal meteors in ref. 9). (This assumption, of course, 
at least partly contradicts the assumption of symmetry upon which was based much of 
the work reported here. ) 

The effect of assumed movement of the j and k planes southward and sunward will 
now be examined. (The consequent modification of observed flux frequencies could be 
computed in a repetition of the entire analysis presented here, but this modification will 
be neglected at this time. ) 

In figure 4, the ordinates at n = 2 for the two sides of the i plane will be assumed 
correct. Any attempt to modify them by assuming that the i plane should be moving 
relative to the Earth would only shift one of them upward and the other downward, de- 
feating the concept of equal flux of kinetic energy in any direction. So to equate the 
fluxes of kinetic energy for the three planes, it will be necessary to move the ordinates 


at n = 2 for the j and k planes downward to match those of the i plane. 

For the j plane, the value of figure 4 is about 204. With = 

1.2107, the corresponding value of is about 170. The v alue of <Pg%orm 

about 14. 8, so that is about 12. 2. The value of = 170 is the second mo- 
ment of V relative to the base value v„^„^ =0. In the same manner as used 

norm norm 

earlier it can be converted to the base = 12. 2 simply by subtracting the square of 

that value, so that 2 ) ~ 18-4. 

For the i plane, the value is about 45. For t he adjusted value of 

<P^orm ^ neither v2^^^^ = 170 nor vgor m(12. 2) = 

18. 4 will do. Instead, some base value Vr, is needed such that v2 ^ = 45, where 

b 'f'g norm(b) ’ 

Vnorm(b) normal velocity component relative to the base V|^. With = 1.2107, 


^normtb) ~ inverse of the process that was used before, the value 

aaT.2) = IS- ^ converted to = 37 by adding (vj, - So, 
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( 12 ) 


~2 _ “2 ~ 
'^norm(b) '^norm(12.2) 


+ (Vh - V, 


norm-' 


or 


37 = 18.4 + (v^ - 12.2)2 


(13) 


and 


= 12.2±4. 3 


(14) 


This result means that, for equating fluxes of kinetic energy between planes j and 

1, the average velocity component relative to the moving plane j must be ap- 

proximately q=4. 3 km/sec rather than the value 12.2. The j plane, then, would need to 
be moving toward the Sun at either 16. 5 or 7. 9 km/sec. The value of cpjv would 

become 5. 2 km/sec. As symmetry of velocity components relative to the j 

plane would now be assumed, the minus values of amount to a reversal of 

flux direction, and the consequent modification of observed influx rates would be ex- 
treme. Hence, the sunward movement of the j plane at 7. 9 km/sec will be chosen for 
use, with the reservation that it should probably be much reduced if the modification of 
observed influx rates had not been neglected. (To avoid reversal of flux, the smaller 
value of V|^ must be used so that the southward movement of the meteoroids may still 
permit them to overtake the j plane.) 

Now let us compare the momentum flux rate for the j plane, with =5.2 

with that for the i plane. From figure 4, for the leading side of plane i, <?^gVj^orm “ 

2. 63 with ^g = 0.311. For the trailing side, 6 with 1.223. The 

average value for the two sides of the i plane, then, is 


. _ 2.63 + 7.6 

^g'^norm « 


5. 1 


(15) 


This result is the same as the adjusted value of 5. 2 for the j plane well within the 
probable accuracy of the determinations. 

A similar analysis for the k plane showed that equalizing the flux of kinetic energy 
with that of the i plane would require an average velocity component normal to the 
moving k plane of 10. 5 or 8. 1 km/sec. For the same reason as before, the value 
8. 1 km/sec was chosen for movement of the k plane in the southward direction with a 
value of = 1-22 km/sec relative to the moving k plane. The value of <PgVj^orm" 
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relative to the moving plane k, then, is 1.22x0.9126 = 1. 11. This rate of southward 
motion of the k plane at 8. 1 km/sec is probably much too high, because reduction of 
the average component of velocity from a value of about 9. 3 derivable from figure 4 to 
the value 1. 22 would mean a very substantial correction to the observed influx rate, that 
is, a substantial reduction of the factor (p . The large reduction in the factor <p would 

o o 

help to bring the kinetic energy flux for the k plane down toward that of the i plane 
without the need of assuming such a large velocity of southward motion for the k plane. 

For the movement toward the Sun at 7. 9 and southward at 8. 1 km/sec the origin for 
the i, j, and k planes, moving at Earth velocity of 30 km/sec, would have to move at 
an angle equal to the arc sine of the expression 



7.9 


8 . 11 ^ 

^ + 8 . 11 


2 


or at an angle of about 21° to the apex of Earth movement. This result does not seem 
unreasonable, particularly in view of the fact it would be considerably reduced when 
taking into account the reduced frequency of observed impacts that would occur normal 
to the j and k planes because of the movement of those planes. 

A rough check relative to existence of the toroidal flux was made. From a study of 
the results of the work of reference 5, it was estimated that the meteoroid sample of 
reference 7 would satisfactorily cover a solid angle defined by the spherical surface be- 
tween two planes, both containing the apex and antapex of Earth movement. One such 
plane would extend through the apex -antapex line, away from the Sun, southward, at a 
dihedral angle of 44° to the ecliptic plane. The other plane would extend through the 
apex -antapex line northward and toward the Sun at a dihedral angle of 44° from the j 
plane . 

Part of the computer program used in this study, as has been described, was re- 
peated under two conditions: (1) with restriction of the sample to meteors arriving from 
the solid angle between two planes extending through the apex -antapex line, both away 
from the Sim, one northward and the other southward, at dihedral angles of 44° from the 
ecliptic, and (2) with restriction of the sample to meteors arriving from the solid angle 
between two planes extending through the apex -antapex line, both northward, one toward 
the Sun and the other away from the Sun, at dihedral angles of 44° from the j plane. 

The first result showed an average component toward the k plane (southward through 
the ecliptic) of 1.65 km/sec. The second result showed an average component of velocity 
toward the j plane (toward the Sun) of 2. 18 km/sec. The component 2. 18 km/sec to- 
ward the Sun is 27. 6 percent of the value of 7. 9 computed to equalize flux of kinetic 
energy. The value 1. 65 km/sec is 20. 3 percent of the value 8. 11 computed as necessary 
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in the southward direction. The components equal to 1. 65 and 2. 18 km/sec would call 
for movement of the origin of the i, j, and k planes at an angle of only about 1/2° to the 
apex of Earth movement, as compared with the earlier computed angle of about 21°. 

The correct values should probably lie between those found by the two methods. 

The calculated values of 7. 9 and 8. 11 for movement of the j and k planes, and the 
angle of 21° relative to the apex, are too high because of neglect of the effect of move- 
ment of the j and k planes on observed values of <p^. The velocity components of 
2. 18 and 1. 65 km/sec, and the angle of about 1/2°, found with the computer run, tended 
to be too low because it was necessary to neglect the directions from which the greatest 
components of velocity in the j and k directions would come, but tended to be too high 
because no correction was made for changes in the values of cp . 

The effect of the movement of the planes on the values of <pg could readily be cal- 
culated in a later analysis. The frequency determined in the present analysis for a 
meteor having a given heliocentric velocity vector would simply be multiplied by the ratio 
of its velocity relative to the moving origin of the i, j, and k planes to the geocentric 
velocity. 

For an ideal analysis it would be desirable to have a repetition of the work reported 
in reference 7 at three additional locations of the cameras: (1) within the arctic circle, 
(2) within the south temperate zone, and (3) within the antarctic circle. The data from 
the additional positions on the Earth's surface would provide coverage of the entire 4ir 
steradians of solid angle surrounding the Earth, except a small solid angle centered 
exactly in the direction of the Sun. Considering the foregoing discussion, it should not 
be surprising if such an ideal analysis with the additional data available would show 
exactly equal flux of kinetic energy across any plane with origin moving at an angle of 
only a few degrees to the apex of Earth movement. 

With the demonstrated existence of such a condition in Earth orbit about the Sun, a 
similar condition could reasonably be expected elsewhere within the solar system, also. 
The sample of reference 7, with similar samples from the other three suggested regions, 
would provide samples for a family of heliocentric vector velocities coming from any lo- 
cation within the solar system. The corresponding velocities and directions at such lo- 
cations are easily computed. As a bimodal log-normal equation has been shown to fit 
the distribution of any velocity component satisfactorily in Earth's orbit, such an equa- 
tion might reasonably be assumed to apply elsewhere also. The maximum possible ve- 
locity and a minimum possible velocity (to avoid passage into the Sun's atmosphere) can 
be computed for any direction of motion at any position within the solar system. 

All the known parameters and conditions described would seem to be sufficient to 
construct a reasonable tentative distribution of velocity and mass and a flux rate for all 
locations. A check would have to be made, of course, as to whether the meteor sample 
included enough asteroidal meteors to indicate that the same findings applied to them. 


21 



DISCUSSION OF RESULTS 


All the work reported here and in references 1 to 5 has culminated in the velocity 

histograms shown in figures 1 to 3; equation (3) with the values of constants shown in 

table I governing distribution of normal components of impact velocity; the concentration 

factors shown in table II; the values of v^ „ shown in figures 4 to 6; and the armor 

norm 

thicknesses shown in table m. The preliminary work reported in references 1 to 5 seems 
to provide a firm foundation for those results. 

The computational results in table HI show a marked advantage of a particular 
orientation only in the case of the trailing side of the i oriented plane. The advantage 
in that case, however, is substantial and can readily be used. 

For example, a vulnerable area of 2000 square feet (185. 8m) may be assumed, 
composed of two parallel plane sheets of armor. Then, if the sheets have the i orien- 
tation with a 99. 5-percent probability of no destructive impact for each sheet (a 99 per- 
cent probability for the combination) the combined thickness for leading and trailing 
sheets, as shown in table HI would be 0. 808 centimeter (0. 318 in.) an average of 0.408 
centimeter (0. 159 in. ). With the j or k orientation, having symmetrical meteoroid 
influx, the combined thickness would be 1. 18 centimeter (0.464 in.) or 1.00 centimeter 
(0. 394 in.), that is, 46 percent or 24 percent greater than with the i orientation. 

If the vulnerable area of 2000 square feet (185. 8 m ) is assumed to be entirely in 
the form of cylinders of uniform wall thickness, with no consideration for shadowing of 
one surface by another, then the necessary thicknesses for the i, j, and k orientations 
are 0. 547, 0.481, and 0. 521 centimeter (0.215, 0. 190, and 0.205 in.). It is seen that 
the i orientation requires 13 percent and the k orientation 8 percent greater thickness 
than the j orientation. 

If cylinders of nonuniform wall thickness are assumed, an approximation of the 
necessary thickness at four points on the cylinder cross section may be obtained by use 
of the results for planes in table m. Those four points are the points of tangency of the 
cylindrical cross section with four planes having i, j, or k orientation as the case may 
be. Thus, an i tube would have two j oriented and two k oriented tangent plane sur- 
faces as illustrated in sketch (b). There would be approximately an area of 1000 square 
feet (92. 9 m ) in each orientation, with a 99. 5 percent probability of no destructive 
impact for each (a combined probability of 99 percent). From the tabdlated results in 
table m, the average wall thickness would then be 0. 545 centimeter (0.2145 in.), with 
0. 50 centimeter required for the tangent k planes and 0. 59 for the tangent j planes. 
Similar treatment for a j or k oriented tube would involve 500 square feet (46. 5m) 
and a 99. 75 percent probability for each of the i planes (leading and trailing) with 
1000 square feet (92. 9 m ) and 99. 5 percent probability for the symmetrical k or 
j planes. The consequent average thicknesses for the j and k oriented tubes would 
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Cross section of i oriented cylinder with exaggerated variation of wall thickness 

(b) 


be 0. 452 and 0. 497 centimeter (0. 178 and 0. 196 in.). These results indicate an ad- 
vantage of 20. 5 and 10 percent for the j orientation as compared with the i or k 
and an advantage of about 7 percent for the j oriented tube of variable thickness over 
the j oriented tube of uniform thickness. 

A further advantage could be obtained by suitable orientation of panels of tubes. 
With use of the most favorable j orientation of tubes of nonuniform wall thickness, a 
large number of such tubes could be placed side by side with all their centerlines lying 
within a single plane, whether i oriented or k oriented. With i orientation of the 
plane of the panel, as shown in sketch (c) below, tubes would, at least to some extent. 


Enclosing lines only to 



23 





TABLE IV. - DATA USED IN CONSTRUCTION OF HISTOGRAMS FOR NORMAL COMPONENT OF IMPACT VELOCITY 

ON PLANES, CYLINDERS, AND SPHERES 


I Velocity, 
km/sec 


Percentage of all sporadic meteors with impact velocity between I and I + 1 


(I) 

Full velocity on 

1 Normal component only 


spherical or 
randomly 

Plane surface whose normal always points 

1 

Cylindrical surface whose axis always 

Spherical 


! oriented plane 


in direction stated 


points in direction stated 

surface 


surface 

1 Direction of movement of I 

Toward Sun 

Perpendicular to 

Direction of 

Toward Sun 

Perpendicular to 




Earth in its orbit 



ecliptic plane 

movement 


ecliptic plane 




Leading side | 

Trailing side 



of Earth in 
its orbit 




0 

0 

6.5732 

1.5325 

0.2150 

1.5865 

1.0611 

1.6925 

1.3149 

1.2706 

1 

.2704 

8. 1473 

3.8174 

1. 1860 

' 4.6238 

3. 1646 

4.9682 

3.5136 

3.6291 

2 

1.0296 

13.4415 

5.9704 

1.9630 

5.7138 

4. 3097 

7.0795 

5.2125 

5.2710 

3 

.5130 

17.2790 

8. 1371 

5.6162 

4.8475 

5.6099 

9.0543 

7.4423 

6.8022 

4 

2.9600 

7.5926 

17. 1491 

4. 8446 

8. 7479 

7.2160 

10.9618 

7.5733 

7.7142 

5 

2.0318 

7.0298 

12.5949 

6.6070 

11.5907 

6.7885 

11.4111 

7.7189 

8.2930 

6 

2.7082 

13.6890 

16.4395 

5.0747 

8.9228 

7.3753 

10.9419 

8. 5386 

9.0178 

7 

9.0513 

4. 1735 

18. 1031 

8. 3602 

8.4837 

7.4244 

9.7357 

8.2444 

8.6125 

8 

6.2843 

1.9316 

10.6094 

7. 8416 

5. 1294 

7. 5947 

7.0078 

8. 7261 

7.6138 

9 

6. 5774 

1. 1193 

4.3381 

7.1637 

8. 1732 

7. 1482 

5.7238 

6.5540 

7.0513 

10 

7. 0204 

1. 7880 

.2665 

5.0596 

4.6995 

5. 7444 

4. 3772 

5.8750 

6. 3905 

11 

9.4131 

1.9782 

.9409 

8.2002 

5. 1469 

6.0317 

3. 8775 

5.3161 

5.2535 

12 

7.7647 

.9120 

.0260 

3.6104 

3. 7602 

5.2222 

2.5981 

3.9572 

4.3015 

13 

7. 1902 

.8027 

0 


5. 7946 

5.0821 

3. 9251 

2. 1453 

3.2478 

3.4277 

14 

6,4739 

1.4562 



2.3588 

1.6657 

3. 1523 

1.6217 

2.6220 

2.6823 

15 

4. 1989 

2.6201 



2.2329 

1.4469 

2. 5991 

.9933 

1.9926 

2. 1309 

16 

3.6807 

.4536 



3. 8782 

.5819 

2.2138 

.6993 

1.9178 

1.7558 

17 

3. 1782 

1.2977 


f 

2.2415 

1.2330 

2.0230 

.7820 

1.5263 

1.4960 

18 

3. 1263 

1. 1581 

.0054 

3.0806 

1.2587 

1. 8489 

.5485 

1.4625 

1.2256 

19 

2.3844 

.7511 

0 


2.2110 

.7287 

1.5393 

.4436 

1. 1954 

1.0076 

20 

1.8870 

.4091 



2. 1668 

.6004 

1.3573 

.4071 

.9025 

.8334 

21 

1. 8402 

.2850 



.9519 

. 8871 

1.0530 

.3697 

.6816 

.6733 

22 

1.4215 

.4104 



.9631 

.9675 

.9019 

.3890 

.6475 

.5754 

23 

1.0350 

.0525 



1.5230 

.6456 

.8333 

.2595 

.6483 

.4832 

24 

1. 3859 

. 1130 



1,5475 

.6217 

.7982 

.2337 

.5631 

.3978 

25 

1.1487 

. 1198 



.9606 

.2853 

.5262 

.1521 

.3980 

.3115 

26 

.7544 

.1154 



1.0116 

.2432 

.4484 

. 1423 

. 3610 

.2514 

27 

.9123 

.0471 

1 


.8020 

.2475 

. 3650 

. 1392 

.2803 

.2060 

28 

. 3758 

.0981 

.0565 

.3746 

.2670 

.2967 

. 1305 

.2049 

. 1626 

29 

.3497 

.3090 

0 


.5786 

. 3673 

.2655 

. 1204 

.1626 

. 1433 







shield each other on the sides tangent with the k plane. (Lines enclosing panels of tubes 
in the sketch (c) are intended only to clarify panel orientation. ) With the panel of tubes 
lying in the k plane (also shown in the sketch) , the tubes would shield each other on the 
sides tangent to the i plane. As the two unshielded k planes require 24 percent greater 
armor thickness than the average of the two i planes, obviously the shielding of the k 
oriented surfaces by the i orientation of the panel of tubes is more desirable than the 
shielding of the i oriented surfaces by the k orientation of the panel. 

If a damage criterion is used that includes as factors and m^ with e 

smaller than in equation (4), with X larger, or both, the value 3. 332 for the exponent of 
'^norm equation (6) is reduced accordingly. The effects of orientation are conse- 
quently considerably altered. Obviously a lower value o f the exponent requires thinner 
armor than a higher value. But, in the function the higher values of n place 

greater weight on the highest values of Thus, with n approaching infinity the 

value of the function (vnorm)^^^ approaches the highest value of for any meteor. 

With n = 1, all weighting in favor of higher values of is gone and we have a 

simple arithmetic average. As n approaches zero, we approach a geometric average 
which is weighted somewhat in favor of lower values of 

From inspection of figure 4, the advantage of the trailing side of the i oriented 
plane is progressively lost with diminishing n until, with values of n less than about 
2.0, it becomes even worse than the leading side of that plane. The leading side of the 
i plane becomes no worse than either side of the j plane when n is about 3. 8, and it 
becomes as good as the k plane with n of about 3.0. From figure 5 the differences in 
effect of the i, j, and k orientations for a tube of uniform thickness are increased as 
n is reduced from the value 3. 332. As n is increased, orientation loses effect when 
n is about 4. 25. With higher values yet, the differences of effect of the three orienta- 
tions reverse, so that the i orientation becomes best. 

An intuitive check on consistency of results has been made by determining that the 
computational result for critical thickness of a sphere as shown in table HI is, within 
1 percent, the same as the average for planes and within 2^ percent the same as the 
average for cylinders with the three principal orientations. A more definite check on 
the result for a sphere is the fact that it is, within 2 percent, the same as determined 
with use of a different method in reference 4. In that case, a histogram was constructed 
with use of the data shown in column 2 of table IV, for gross impact velocities. An 
equation of the same form as equation (3) was fitted to that histogram and that equation 

was integrated to p rovide values of v”. Then, with cp^ equal to unity, and with use of 

an expression for in terms of v^, the value of t^^ was determined with an 

equation equivalent to (6) . 

The analysis here has dealt entirely with conditions far enough from Earth that 
Earth’s gravity has negligible effect on meteoroid velocity. The condition at the upper 
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limit of Earth’s atmosphere was also treated, for a sphere, in reference 4. The neces- 
sary armor thickness for the sphere near Earth proved to be about 19 percent greater 
than far from Earth. 

The results for the near -Earth condition for a sphere were simply derived in ref- 
erence 4 from the far -from -Earth results. But the treatment of effects of orientation of 
planes and cylinders for near-Earth conditions would be much more complex and could 
not be included in the present analysis. The author believes that the relative advantages 
of particular orientations would remain of the same nature as for far -from -Earth con- 
ditions, but that they would be slightly reduced. This opinion is based on the fact that 
substantial changes of direction of motion of meteoroids by Earth's gravitation exists 
only when values of Vq are too small for a meteoroid to represent a substantial part of 
the total hazard. 

Although the results of this study are directly applicable only to space missions in or 
near Earth's orbit about the Sun, it is believed the results can be extended by the methods 
outlined to deduce a valuable tentative model at substantial distances from that orbit. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, Jime 22, 1970, 

131-05. 
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APPENDIX A 


SYMBOLS 


A 

^eff 


a 

b 


Ci 


c 


®e(s) 

®k(E) 


e 


2 

vulnerable area exposed to meteoroid hazard, m 

projection of vulnerable area exposed to meteoroid hazard upon a plane nor- 
mal to meteoroid path, unspecified units 

semi-major axis of elliptical orbit of meteoroid, or constant as used in stand- 
ard eq. (2) for hyperbolic orbit, AU 

semi-minor axis of elliptical orbit of meteoroid, or constant as used in stand- 
ard eq. (2) for hyperbolic orbit, AU 

constant in eq. (3), dimensionless 

an orbital element for an elliptical, parabolic, or hyperbolic orbit, as illus- 
trated in fig. 7(a), AU 

energy of escape from Sun at positions in Earth's orbit, ft -lb (1. 356 J) 

kinetic energy associated with orbital motion of a mass corresponding to one 
unit weight in Earth's orbit, ft -lb (1. 356 J) 

eccentricity of conic section (orbit of meteoroid), dimensionless 


Fgs^p) gravitational force exerted by Sun on mass corresponding to unit weight on 
Earth's surface, but at distance P from Sun, lb (4.448 N) 

F(Zd)„-, a statistical function developed in reference 1 involving zenith angle of meteor, 

XV ciV 

its position within camera field of view, and azimuth of its path, and their 
effects on density of photographic trace, unspecified units 

f( ) (with any argument within parentheses) statistical frequency of occurrence of 

given value of argument (for , approximately fraction of all impacts 

for which has given value ±0. 5 km/sec) 

2 

g acceleration due to Earth's gravity at Earth's surface, ft sec (0. 3048 m 

sec“^) 

i direction of the apex of Earth's movement, also inclination of orbit of meteor- 

oid to ecliptic plane as used in ref. 7 (equal to y or supplement of y as 
used here), deg 

r unit velocity vector in the i direction, km/sec 
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^cyl 

Vl 

j 

J* 

k 

k 

m 

P 

Pe 

q 

q’ 

^cyl 

^sph 

®cyl 

^sph 

T 



for meteoroid that passes through axis of i, j, or k oriented cylinder, angle of 
incidence of meteoroid’s path relative to axis of cylinder, deg 

angle of incidence of meteor path relative to plane with i, j, or k orientation, 
deg 

direction within ecliptic plane and approximately away from Sun but exactly per- 
pendicular to apex of Earth motion 

unit velocity vector in j direction, km/sec 

direction perpendicular to ecliptic plane, southward 

unit velocity vector in k direction, km/sec 

mass of meteoroid, g 

average number of destructive impacts by meteoroids per mission or semi-latus 
rectum of meteor orbit, AU 

sec per year 

perihelion distance of meteoroid orbit, AU 

aphelion distance of meteoroid orbit, AU 

radius of i, j, or k oriented cylinder, unspecified units 

radius of Earth's orbit, ft or AU 

radius of sphere, unspecified units 

cylindrical surface area for an i, j, or k oriented cylinder, unspecified units 
surface area of sphere, unspecified units 
duration of space mission, sec 

critical value of thickness of armor necessary to provide a given probability of 
no destructive impact by a meteoroid, cm 

velocity of Earth in its orbit, km/sec, except where shown as ft/sec 

velocity of escape from Sun at positions in Earth’s orbit, km/sec except where 
shown as ft/sec 

scalar value of velocity of meteoroid relative to moving Earth, corrected for 
acceleration by Earth gravity, as given in ref. 7, km/sec 

vector velocity of meteoroid relative to moving Earth, corrected for acceleration 
by Earth's gravity, km/sec 

scalar value of velocity of meteoric particle relative to Sim, corrected for ac- 
celeration by Earth’s gravity, as given in ref. 7, km/sec 
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V 


norm 



6 . 

1 

€ 

h 

e 


vector velocity of meteoroid relative to Sun with scalar value Vjj, km/sec 

component of the vector velocity of a meteoroid ^ having scalar value I and 
normal to surface of cylinder with i, j, or k orientation, km/sec 

(with subscripts i, j, or k) component of vector ^ in i, j, or k direction, 
km/sec 

(with subscripts ir, jr, or kr) component in i, j, or k direction of vector 
velocity of meteoroid relative to moving Earth v^, km/sec 

(with subscripts ij, jk, or M) component lying entirely within k, i, or j 
oriented plane, respectively, of vector velocity of meteoroid relative to moving 
Earth km/sec 

component of in i direction under the usually false assumption that meteor- 
oid orbit is within ecliptic plane, km/sec 

scalar value of component of impact velocity of a meteoroid normal to impacted 
surface, km/sec 

scalar velocity of meteoric particle relative to Earth's atmosphere before de- 
celeration by atmosphere, km/sec 

one of rectangular coordinates in standard equations for conic sections, AU 
one of rectangular coordinates in standard equations for conic sections, AU 
angle of meteor path through atmosphere relative to zenith, deg 
right ascension of true radiant of meteor (ref. 7), deg 

inclination of orbit of meteoroid to ecliptic plane (arbitrarily measured with 
value equal to or less than 90) , deg 

declination of true radiant of meteor (see ref. 7), deg 

value of angle at which normal component of impact velocity of a particular 
meteoroid upon an i, j, or k oriented cylinder or upon surface of a sphere 
is equal to I, deg 

offset of mode in i^^ term of log-normal distribution eq. (3), km/sec 
exponent of in damage criterion 

percent elongation of armor material 

angle between apex of Earth motion and line within ecliptic plane that is per- 
pendicular to major axis of orbit of meteoroid under usually incorrect assump- 
tion that meteoroid’s orbit is within ecliptic plane (see fig. 7(a)), deg 
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\ exponent of m in damage criterion or angle between vector velocity of meteor- 

oid and apex of Earth motion, imder usually incorrect assmnption that 
meteoroid orbit is within ecliptic plane, deg 

length of an i, j, or k oriented cylinder, unspecified units 

logarithmic modal value in term of distribution eq. (3), natural logarithm of 
velocity in km/sec 

P distance from Sim, ft (m) 

Q 

P density of particle impacting upon armor, g/ cm 

Jr 

O 

density of armor material, g/cm 

CTj logarithmic standard deviation in a term of distribution eq. (3), natural log- 

arithm of value of velocity in km/sec 

cp angle between vector velocity of meteoroid and line within ecliptic plane 

that is perpendicular to major axis of orbit of meteoroid under usually false 
assumption that meteoroid’s orbit is within ecliptic plane (see fig. 7(a)), deg 

(pcyi()i (with i, j, or k within parentheses and with integral value of I) 

^cyl(j)’ ‘^cyl(k) fraction of cases in which given meteoroid should be 

expected to strike cylinder with i, j, or k orientation with normal component 
of impact velocity between 1 and I + 1 km /sec 

‘^cyl(m)n " ‘Pcyl(m) significance of m as i, j, or k^ 






^g(qu) 


cPl 


concentration factor, frequency of meteoroid impacts upon structure under 
specified condition relative to frequency of impacts upon same structure ran- 
domly oriented 


summation of values of for particular significance of m and particular 

significance of n 


same as cpg, but applying only to meteoroids for which v- • j was negative 
and v^ • ^ was positive 


fraction of <Pcyi(i)’ ‘^cyl(j)’ ‘^cyl(k) which normal component of impact 
velocity of meteoroid upon cylinder is between I and I + 1 or fraction of 
(p^ for which normal component of impact velocity on sphere is between I 
and I + 1 
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^m(n) 

*^111 (r) 
^or 

^ or(av) 
‘^sph(I) 

‘^w(r) 


(m signifying pi or cyl and n signifying i, j, or k) a modification of the 
weighting factor (p^ for use with a particular meteor, for a plane or cylinder 
with i, j, or k orientation, to take into account reduced effective vulnerable 
area because angle of incidence is not the most unfavorable. 

(m signifying cyl or pi) same as (p^ but referred to randomly oriented 
cylinder or plane 

probability of impact of given meteoroid upon given structure (plane or cylinder), 
with i, j, or k orientation, divided by probability of impact upon same struc- 
ture with basic (worst) orientation relative to meteoroid path 

average value of <p^^ for given structure with given orientation 

weighting factor for correction of spacewise bias in photography of meteors in- 
verse value of which was developed and computed in work reported in ref. 5 

that part of (p^ applying to impact upon sphere with normal component of 
impact velocity between I and I + 1 km/sec 

weighting factor developed in refs. 1 to 5, to provide simultaneous correction 
for mass, velocity, and spacewise biases in photography of meteors 

same as ip^ but referred to randomly oriented area rather than area always 
normal to meteoroid path 
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APPENDIX B 


DATA USED FOR EACH METEOR 

Redundancy exists in the orbital parameters of meteors given in reference 7. Some 
of those parameters are (1) the angle 6 of the true radiant, (2) the angle a of the true 
radiant, (3) the velocity of the meteoroid relative to the Sun, Vjj, corrected for the ac- 
celeration by Earth gravitation, (4) the semi -major axis a, (5) the eccentricity e, 

(6) the perihelion distance q, (7) the aphelion distance q’, and (8) the inclination of the 
meteoroid's orbit to the ecliptic plane y (designated i in ref. 7). 

Because of symmetry assumed in this analysis, unnecessary and useless compli- 
cations would be Involved in making use of the values of 6 and a, or any of several 
other parameters given in reference 7 and not listed here. Of the parameters listed, 
three only are needed, including the orbital inclination y. Values of a, q, and y were 
used in determination of angles defining the direction of movement of the meteoric par- 
ticle relative to the Earth's orbital plane and the apex of the Earth's movement within 
that plane. The value of Vjj as given in reference 7 was used directly as the velocity 
of the particle. 

Arrangements were made within the computer program to use the values of a and 
q to compute the eccentricity e, the velocity Vjj, and the aphelion q'. These values 
were computed by well known methods that will not be described here and were com- 
pared for all meteors used with the values given in reference 7. No substantial dis- 
agreements were found. 
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APPENDIX C 


BASIC CONSTANTS USED IN RESOLUTION OF VELOCITIES OF METEORS 

Certain basic constants will be explained here, which were needed in the process of 
solving for velocity components, beginning with given orbital parameters. For conven- 
ience, the customary U.S. system of units was used in this appendix and the results to 
be used in the analysis were converted to SI imits after the determinations were com- 
pleted. The obvious constants will not be derived. 

The radius of Earth's orbit was taken as 

rg = 4. 9106X10^^ ft (14. 9714 km) (Cl) 

The length of year was taken as 

Pg = 3. 1558X10 sec (365. 256 days) (C2) 

The Earth orbital velocity was taken as 

Vg = = 9. 7770x10^ ft/sec (32. 069 km/sec) (C3) 

Pe 

The acceleration due to Earth's gravity at surface, at the camera sites in New Mexico, 
was 


g = 32. 2 ft/sec^ (9. 815 m/sec^) 


(C4) 


The kinetic energy of a mass corresponding to one unit weight in Earth’s orbit is 


Ej^^gj = ± v| = 1. 4843X10^ ft -lb (2. 011x10® J) 


(C5) 


The velocity of escape from the Svin at a position in Earth’s orbit was taken as 


Ve(s) = \^Vg = 1. 3827x10® ft/sec (42. 14 km/sec) 


(C6) 
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Equation (C6) is derived by integrating the Sun's gravitational force from a distance 
of one astronomical unit to infinity. The escape kinetic energy must be 


E 


e(s) 


-r 


^E 


’gs(p) 


dp 


(C7) 


where is the force exerted by the Sun on a mass corresponding to one unit weight 


at any distance p. An expression for F 


gs(P) 


may be derived from the equality of cen- 


trifugal force and gravitational force of a particle moving in the Earth's orbit as follows: 


'gs(rg) 


1 '"e 


(C8) 


g r 


E 


and as the Sun's attraction varies inversely as the square of the distance, 


' gs(P) 


,,2 ,2 

_ 1 ^E / e\ 

g r g V P / 


(C9) 


From equations (C7) and (C9), 


r - 1 v2 r 
^e(s) - - '^E^'e 


/ 


p‘2dp = i v| 


(CIO) 


■E 


On comparison of equations (C5) and (CIO) we see that E^^g^ is twice So equa- 

tion (C6) follows because velocities must vary as the square roots of kinetic energy. 



APPENDIX D 


THREE ORTHOGONAL COMPONENTS OF METEOROID VELOCITY 

RELATIVE TO EARTH 


Figure 7 illustrates the method of resolution of the heliocentric velocity of a meteor- 
oid Vjj into three orthogonal components. As an intermediate step, the orbit of the 
particle was treated as being within the plane of the ecliptic (fig. 7(a)). Although two 
points of intersection always exist, only one need be considered because of symmetry 
in the use that is to be made of the results. The upper intersection in figure 7(a) was 
chosen for consideration. 

Also because of symmetry, no discrimination was necessary between positive and 
negative j and k, directions. In figure 7(a), for convenience, the j direction was 
taken as the reverse of the actual for the meteors in direct orbits that were treated. 
Moreover, at this stage of the analysis, all meteors were treated as having direct 
motion, that is, angular velocity about the Sun in the same sense as for Earth. Dis- 
crimination between direct and retrograde movement in the orbit was left till later. 

The heliocentric vector velocity of a meteoroid will be designated as Then, 

of course. 



(Dl) 


In the intermediate step, this vector velocity ^ was resolved into two components ’vT 
— ^ ^ 3 

and Vj. The component Vj was the final component wanted in the j direction, but the 

component vJ was simply the vector velocity ^ minus 'vT. It was the component ve- 
locity in the i direction that would have existed if the particle’s orbit had actually been 
in the ecliptic plane as shown in figure 7(a). 

Obviously, in the intermediate step. 


Vg sin xT 


and 



cos Xi 


(D2) 


(D3) 


where i and j are vuiit vectors in the i and j directions, and X is the angle be- 
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tween the tangents to the meteoroid path and Earth’s path in their orbits at the time of 
impact (fig. 7(a)). 

Methods used for finding the value of the angle X will now be described for ellip- 
tical, parabolic, and hyperbolic orbits. Sources of various equations will not be stated 
in cases where they are well known and are obtainable in most textbooks on analytic 
geometry. 


Elliptical Orbits - Angle Between Paths of Meteoroid and Earth 

If c represents distance from the center of an elliptical orbit to the applicable 
focus shown in figure 7(a), then, obviously, 

c = a - q (D4) 


Then the semi-minor axis of the orbit is 

b = ^a^ - c^ 


(D5) 


and the eccentricity is 


e = (D6) 

a 

Now the angle 9 with vertex at the Sun in figure 7(a) may be found with use of the 
standard polar equation of an ellipse about its right-hand focus as a center. 


X p, — — 

1 + e cos 9 

and the equation for the semi-latus rectum 


(D7) 


P 



(D8) 


(r g is unity in eq. (D7) and elsewhere with distances expressed in astronomical units) . 
By combining equations (D7) and (D8) and rearranging. 
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- r„a 

cos 0 = ^ (D9) 


Next, the slope of the tangent to the particle orbit at the point of intersection of the 
orbits was determined. By implicit differentiation of the equation of the ellipse 


the slope is foimd as 



dx 



or 


dx 


(c + r g cos 0) 
a^r j, sin 0 


The angle cp as shown in figure 7(a), then, is 


(p 



IL 

2 


At this point, the angle X (fig. 7(a)) is determined as 


(DIO) 


(Dll) 


(D12) 


(D13) 


X = 0 - (p (D14) 

Parabolic Orbits - Angle Between Paths of Meteoroid and Earth 

The equation for a leftward extending parabola like that shown in figure 7(a) is 

= 4 cx (D15) 
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in which c must be negative. As may be seen from figure 7(a) , 

2 / ^2 2 
rg = (c - x) + y 

From equations (D15) and (D16) 

rg = - 2cx + + 4cx 

or 


-r j. = c + X 


and 


X = -rg - c 


As, for the parabola, 


c = -q 


an expression for 6 may be written by reference to the figure. 


a q + X 
cos 6 = -1 


or, from equations (D19), (D20), and (D21), 


cos 9 = - 1 


The angles 9 and X may now be foimd as follows: From equations (D 15), 
and (D20), 

y2 = -4q(-rg + q) 


or 


(D16) 

(D17) 

(D 18 ) 

(D19) 

(D20) 

(D21) 

(D22) 

(D19), 

(D23) 
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By implicit differentiation of equation (D15), with substitution according to equation (D20) , 


(D25) 

dx y 


From equations (D24) and (D25), 


dx 



(D26) 


The angles cp and X may now be found with equations (D13) and (D14) as before. 


Hyperbolic Orbits - Angle Between Paths of Meteoroid and Earth 

Equation (2) was used for the hyperbolic orbits. The variable a as used in the 
denominator in the first term of the equation was made equal to the absolute value of a 
as given in reference 7. For the leftward branch of the hyperbola as shown in fig- 
ure 7(a), the value of x determinable with equation (2) is always negative. 

As shown in figure 7(a), for the hyperbola, 

c = q + a (D27) 


Hence, with the well known relation. 


c= 


(D28) 


the value of b may be expressed as 


b = |/(2a + q)q (D29) 

With use of the polar equation for the leftward branch of a hyperbola about the left- 
ward focus, the condition in figure 7(a), 



In equation (D30) e and p are expressed by equations (D6) and (D8) for the hyperbola 
as they were for the ellipse. Then, from equations (D6), (D8), and (D30), 


cos $ = 






By implicit differentiation of equation (2), 


dy _ b^x 


By inspection of figure 7(a), 


X = -c + rg cos 6 


and 


y = rg sin 9 


(D31) 


(D32) 


(D33) 


(D34) 


From equations (D32) to (D34), 


dx 


b^(r 


E 


a^r 


cos 0 - c) 


E sin 9 


(D35) 


After determination of dy/dx with equations (D32) to (D34), the angles cp and X 
may be found with use of equations (D13) and (D14), as with the elliptical and parabolic 
orbits. 


Parts of Heliocentric Velocity Resolution Common to All Orbital Types 


Because of the simplifying approximation of Earth's orbit as a circle, application 



of orbital equations was not possible for meteors having a value of q greater than 1 as- 
tronomical unit or values of q’ less than such a value. For all such meteors, the value 
of q or q' was taken as unity, and impact with Earth's atmosphere was assumed to 
have occurred exactly at the point of tangency of the orbits of the particle and the Earth, 
under the presumption used in this intermediate step that the orbit of the meteoroid lay 
entirely within the ecliptic plane. Hence, for such cases, as well as for cases in which 
the value of q or q' was given in reference 7 as exactly 1.00, no discrimination was 
needed in accordance with the preceding sections between elliptic, parabolic, and hyper- 
bolic orbits. Instead of equations (D2) and (D3), the following were used 

■^=VjjT (D36) 


and 


TT= 0 (D37) 

After determination of the values of v^ and "v^ with use of equations (D2) and (D3) 
or equations (D36) and (D37), the second and final stage in determining the three prin- 
cipal components of v^ was performed as illustrated in figure 7(b). The relation of 
this view to figure 7(a) is shown in figure 7(a); that is, figure 7(b) shows the orbits as 
they would be seen by a viewer looking in the direction from Sun to Earth. In figure 7(b), 
however, the orbit of the meteor particle has been rotated about the line that extends 
from Sun to Earth through the angle y. This angle is the same as i, given in refer- 
ence 7, except that when the particle was in retrograde motion and the value of i was 
consequently shown as greater than 90°, the angle y was made equal to the supplement 
of i. 

As is obvious by inspection of figure 7(b), 

1^1 sinyk (D38) 

where k is a unit vector in the k direction (normal to both i and j introduced earlier), 
and 




cos yi 


(D39) 


From equations (D3) and (D38), 
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(D40) 



cos X sin yk 


and from equations (D3) and (D39), 


Vj = Vjj cos X cos yi (D41) 

Equation (D2) or (D37) and equations (D40) and (D41) provide the desired three principal 
components of the velocity for particles in direct motion. Equations (D2), (D37), 
and (D40) are also valid for particles in retrograde motion because of symmetry. For 
retrograde motion, equation (D41) is valid with use of -i in place of T. 


Velocity Components Relative to Moving Earth 


From the values of ”v^, 'vT, and v^, and the velocity of Earth in its orbit v„, com- 

1 J -Cj 

ponents of the vector velocity of the meteoric particle relative to the moving Earth 
were next determined. The vector velocity v^ always had approximately the same 
scalar magnitude as the velocity v^ given in reference 7. The components of v^ were 
determined with the equation 



(D42) 


and with the two equations 



(D43) 



(D44) 
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TABLE V. - COMPUTED DATA FOR THREE SAMPLE METEORS 



Reference or equation 

Meteor number 


number 

6102 

6430 

8761 

Page 

ref. 7 

35 

39 

34 

Line 

ref. 7 

22 

9 

17 

Direct or retrograde 

— 

Direct 

Direct 

Retrograde 

Orbital type 

— 

Elliptical 

Parabolic 

Hyperbolic 

y, deg 

— 

13 

68 

42 

Cos Zj^ 

ref. 7 

0.81 

0.89 

0.95 

v^. km/sec 

ref. 7 

28.7 

43.6 

69.9 

F(Zj^) 

av 

ref. 1 

0.6272 

0.4734 

0. 3863 

a, AU 
q, AU 

ref. 5 
ref. 7 
ref. 7 

0.8581 

2.75 

0.49 

0.6389 

0,4713 

5.87 

0.89 

0. 96 

Vj^, km/sec 

substituted for Vj^, km/sec 

ref. 7 

38. 5 

42.4 
42. 14 

44. 2 




Vq, km/sec 
b (AU) 

ref. 7 

eqs. (D5) or (D29) 

26. 7 
1.567 

42. 1 

68. 9 
3. 352 


c (AU) 

eqs. (D4), (D20), or (D28) 

2. 260 

-0.960 

6. 761 


eq. (D6) 

0. 822 


1. 152 



Cos f) 

eqs. (D9). (D22), or (D31) 

-0. 1209 

0.9200 

0. 7945 

0, deg 

— 

96.94 

23.07 

37. 39 

dy/dx 

eqs. (D12), (D26), or (D35) 

-0.6995 

-4. 8990 

-3. 2046 

CO, deg 

eq. (D13) 

55.03 

11. 54 

17. 33 

A, deg 

eq. (D14) 

41. 92 

11. 54 

20.06 

(■^ 1 , km/sec 

eq. (D2) 

25.72 

8.427 

15. 163 

|v! 1 , km/sec 

eq. (D3) 

28.65 

41. 30 

41. 51 

1 ^ 1 , km/sec 

eq. (D40) 

6.443 

38. 28 

27, 78 

■ T, km/sec 

eq. (D41) 

27.92 

15.47 

-30. 84 

« i, km/sec 

eq. (D42) 

-1. 884 

-14. 33 

-60.65 

. km/sec 

eq. (D43) 

25. 72 

8.428 

15. 163 

|\^L km/sec 

eq. (D44) 

6.445 

38. 28 

27. 78 

, km/sec 

eq. (E3) 

25. 79 

16.626 

62. 52 

|i^!, km/sec 

eq. (E4) 

26. 51 

39. 20 

31.65 

km/sec 

eq. (E5) 

6. 714 

40.88 

66.71 


eq. (1) 

1. 997x10"® 

4. 262x10"'^ 

7.931X10"® 

"^cyl(i) 

eq. (E6) 

1. 992x10"® 

4.003X10"'^ 

3.669X10"® 

^cyl(j) 

eq. (E6) 

5.044X10"'^ 

4. 175X10"'^ 

7. 734X10"® 

"^cyl(k) 

eq. (E6) 

1. 937X10“® 

1.698X10"'^ 

7. 248X10"® 

"’pl(i) 

eq. (E2) 

1.415X10“'^ 

1.464x10"'^ 

7.032x10"® 


eq. (E2) 

1. 932X10"® 

8.607x10"® 

1. 758X10"® 

‘^pl(k) 

eq. (E2) 

4. 842x10""^ 

3.910x10"'^ 

3.221x10"® 

^cyl(r) 

eq. (F17) 

1. 568X10"® 

3. 348X10"'^ 

6.229X10"® 


eq. (F7) 

9. 984x10"'^ 

2. 131x10"'^ 

3.966X10"® 

^cyl(i)0 

eqs. (ElO) and (Ell) 

1.417x10"^ : 

1. 303x10'^® 

1.832x10"^^ 

"^"cyl(i)l 

eqs. (ElO) and (Ell) 

4.258X10"^ : 

3. 911X10"^® 

5. 501X10"^^ 

^cyl(i)24 

eqs. (ElO) and (Ell) 

1. 831x10"'^ 

8. 178X10"® 

1-419x10"® 

^cyl(i)25 

eqs. (ElO) and (Ell) 

2. 731x10"'^ 

8. 749X10"® 

1.579X10"® 

^cyl(i)26 

eqs. (ElO) and (Ell) 

3.904x10“'^ 

9. 371x10"® 

1. 778x10“® 

^cyl(i)79 

eqs. (ElO) and (Ell) 

0 

0 

0 



APPENDIX E 


USE OF ORTHOGONAL COMPONENTS OF METEOROID VELOCITIES 
RELATIVE TO EARTH TO OBTAIN DISTRIBUTIONS OF 
NORMAL COMPONENTS OF IMPACT VELOCITY 

The two procedures outlined in the sections entitled Determination of Velocity Com- 
ponents Normal to Surfaces of Structures from Principal Components of Velocity Rela- 
tive to Earth and Normal Velocity Distributions as Normalized Sums of Adjusted Unit 
Counts will be explained in detail here without a formal separation of the two procedures. 
The operation of the computer program in determining the distribution of normal impact 
velocity will be described on the basis of analogous hand tabulations, tables V and VI. 

As the procedure is simplest for impact upon planes, that case will be treated first. 


Distribution of Normal Component of Impact Velocity Upon Planes 

For planes an assumption was made that impacts would be statistically symmetrical 
relative to the two sides of the j and k planes. No such assumption was possible for 
the i plane. For the meteors treated, all impacts were on the same side of the j 
plane, and all on the same side of the k plane. Accordingly, velocity distributions could 
be determined for impacts on one side of each of these planes and could be assumed to 
apply to the other side also. A separate distribution was obtained for each side of the i 
plane. Negative values of v7^ • T indicated impact on the leading side of the i plane; 
positive values the reverse. 

Table VI(a) is a skeleton tabulation for the j plane. It shows entries for three 
specific meteors and for several specific normal components of impact velocity. In 
table V appear identifications of the three sample meteors, the pertinent basic data as 
given in reference 7 or as calculated in the work of references 1 and 5, the results of 
various intermediate values calculated, and final values entered in the tabulation 

of table VI(a). The entire procedure for arriving at the value specifically for 

meteor 6102 will now be explained. (Table V also includes results for cylinders and 
spheres, which will be discussed later.) 

The sources of the first 29 items in table V are obvious, have been clearly stated 
in the table, or have been explained in the main text or in preceding appendixes. Later 
items pertinent to the tabulation for the j plane in table VI(a) are v^, cp^ from equa- 
tion (1), and 
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TABLE VI. - SKELETON TABULATION OF ADJUSTED UNIT COUNTS RELATIVE TO NORMAL 


COMPONENT OF IMPACT VELOCITY OF METEORS 

[Entries for most values of I and all entries for meteors other than 6102, 6430, and 8761 
have been omitted. ] 


(a) Plane whose normal points toward Sun 


Meter serial 

Fully adjusted unit count for normal component of impact velocity 


Total of 

Fully adjusted 

number 

1=0 

1=7 

between I a 
1=8 

nd (I + 1), km 
1=9 

/sec 

1=15 

1=25 

1=79 

columns 
I = 0 to 
I = 79, 

unit count for 
impact on 
randomly 

6102 

0 

0 

0 

0 

0 

1.932x10"® 

0 

1.932X10'® 

oriented plane, 
^^pl(r)^ 
9.984x10"'^ 

6430 

0 

0 

8.607X10'^ 

0 

0 

0 

0 

8. 607x10"® 

2. ISlxlO""^ 

8761 

0 

0 

0 

0 

1. 758x10"® 

0 

0 

1. 758x10"® 

3.966x10"® 

Totals 

9. 853x10"^ 

3. 831x10"^ 

3. 594x10""^ 

3.283x10"'^ 

1.081X10"^ 

4.402X10"® 

0 

^4. 583X10"® 

^3.786X10"® 

Percentage totals 

.2150 

8. 360 

7. 842 

7. 164 

2.233 

.9606 _ 


100.00 



^Total impacts as percentage of statistically expected impacts on a randomly oriented plane surface = 
(4. 583xl0"^)/(3. 786X10'^) x 100 = 121.07 percent. 


(b) Cylinder whose axis points in direction of Earth movement 


Meter serial 
number 

Fully adjusted unit count for normal component of impact velocity 
between I and (I + 1), km /sec 

Total of 
columns 
1=0 to 1=79, 

Fully adjusted 
unit count for 
impact on ran- 
domly oriented 
cylinder, 

^^cyl(r)^ 


1=0 

1=1 

1=24 

1=25 

1=26 

1=79 

6102 

1.417X10"® 

4. 258X10"® 

l-SSlxlO""^ 

2.731X10"'^ 

3.904x10"^ 

0 

1.992x10"® 

1. 568x10"® 

6430 

1. 303x10"^® 

3.911X10"^® 

8. 178X10"® 

8.749X10"® 

9.371x10"® 

0 

4.003X10"^ 

3.348X10"'^ 

8761 

1.832X10"^^ 

5. 501x10"^^ 

1.419X10"® 

1 . 579x10"^ 

1.778X10"® 

0 

3.669X10"® 

6.229X10"® 

Totals 

6.959X10"® 

2.076x10"^ 

5.235X10"® 

3.451x10"® 

2.941X10"® 

0 

^6.559X10"® 

^5.947x10"® 

Percentage totals 

1.061 

3. 165 

.7982 

.5262 

.4484 

Q 

100.00 





^otal impacts as percentage of statistically expected impacts on a randomly oriented cylindrical surface = 
(6. 559x10 ”^)/(5. 947X10"^) X 100 = 110. 30 percent. 
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The adjusted unit count cp^ must be modified for application to planes (and to cyl- 
inders, but not to spheres) to provide for the effect of variation of effective vulnerable 
area with the angle of incidence of the meteoroid. The modified adjusted vmit count 
(ppj(j) is the result for the j plane. The modification is illustrated for the j oriented 
plane (and for later discussion of the j oriented cylinder) in figure 8. 

In figure 8 the vector velocity of the meteoroid relative to the moving Earth is v^, 
which generally is not in the ecliptic plane. The projection of v^ upon the j plane is 
designated as and its projection on the axis normal to the j plane is v^ (eq. (D43)). 
Now a unit of the j plane will project upon a plane oriented normally to v^ as an effec- 
tive vulnerable area equal to the cosine of the angle i^j, which is expressed as 


Vl 


= tan 



(El) 


As the weighting factor (p^ provides for the correct relative frequency of impact for 
angles of incidence always arbitrarily zero, the relative frequency of impact for other 
angles of incidence should be equal to the effective vulnerable area times <p^. The 
modified adjusted unit coimt for the j plane should therefore be equal to the value of 
(p^ multiplied by the cosine of the angle ipj, or 



The velocity components v^, v^j^, and Vj^, scalar magnitudes for which are needed 
in equations (E2), are obviously obtainable with use of the values of v.^, vT^, and v^ 
(eqs. (D42) to (D44)), as 
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'^ij ''^ir "'^jr 

(E3) 

1 

(E4) 

^jk = '^jr + ^kr 


(E5) 

"^ki = '^kr + "^ir 


Now the modified adjusted unit count for the j plane for meteor 6102 was 

1. 932x10”^ as shown in table V. The scalar magnitude of the normal component of im- 
pact velocity, v.* , was 25. 72 km/sec. Hence, the modified adjusted unit count of 

ft J i 

1. 932x10"° should properly be entered in table VI(a) in the column for I = 25. In a 

-8 -8 

similar manner, the values 8.607x10 and 1. 758X10" for meteors 6430 and 

8761 were entered in the proper columns for | v^^ | equal to 8. 427 and for |vT^| equal 
to 15. 163, that is, in the columns for 1=8 and I = 15. 

In the full tabulation, of which table VI(a) is only a skeleton, all meteors that were 
treated involved the same procedure, with entry of a value of fhe column having 

a value of I determined from the value |v^|. In addition, an entry was made in the 
column for ’’Total of columns 1 = 0 to I = 79. ” This entry was always the same as 
seen later, the procedure was somewhat different for cylinders and 
spheres. (The last column of table VI(a) and the footnote may be ignored until values of 
(p^ are discussed in appendix F. ) 

After all meteors meeting the necessary conditions had been so treated, the columns 
for I = 0 to 79 were totaled, and the results were converted to percentages as shown in 
the last line of the table. The percentage results were entered in the appropriate column 
of table IV and were plotted as the circular symbols in figure 1(c). 

The i and k oriented planes were treated in the same manner, in tabulations 
similar to that represented by table VI(a), with use of values of 
values of vj^ • i and |v^| from table V. For the i oriented plane, however, a 
separate tabulation was made for each side. Negative values of v-^ • T called for entry 
of ill ih® appropriate column of the tabulation for the leading side of the plane, 

positive values the opposite. The resulting velocity distributions for the two sides of 
the i plane and for the k plane appear in table IV and are plotted as the circular sym- 
bols in figure 1(a), (b), and (d). 

As has been explained, only meteors reaching Earth from the side of the j plane 
opposite the Sun and from the north side of the k plane have been considered. However, 
any meteoroid arriving from other directions symmetrically relative to either or both of 
those planes would strike the i, j, or k plane with the same normal component of rela- 
tive velocity. Hence, with the assumption of symmetry of meteoroid influx relative to 
the j and k planes, the method of determining distribution of normal components of 
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impact velocity that has been described, when applied to the one quadrant of space that 
is totally covered, should be valid for impact from all directions. This same statement 
will apply to the method used for cylinders. 


Distribution of Normal Component of Impact Velocity Upon Cylinders 


The manner of determination of distributions of normal components of impact ve- 
locity upon cylinders was similar to that for planes. However, the i, j, and k oriented 
cylinders could all be treated in the same manner, without the need for distinction be- 
tween two sides of any of them. The desired result, for each orientation of a cylinder, 
applies to all parts of the surface collectively. 

A skeleton tabulation for an i oriented cylinder appears in table VI(b). It is just 
like table VI(a) and is made up from data given in table V in a manner that will now be 
described. 

In table V the items to and including have been explained. The items following 
cp^ that are pertinent to the distribution of normal component of impact velocity upon 
the i cylinder are ^c>cyl(i)’ *i^cyl(i)n equal to 0, 1, 24, 25, 26, and 79. 

The modified adjusted unit count <;c’cyl(i) the same purpose for the i cylinder 

as serves for the j plane. This analogous modification of the adjusted unit count 

for the i oriented cylinder may be referred to the worst case, in which Vq would lie 
within the i plane. With reference to such a case, as may be seen in figure 8, the 
effective vulnerable area presented by the cylinder would vary as the cosine of the angle 
icyi> or as the sine of ipj given for the case of the j plane by equation (El). Hence, 


‘^cyl(i) 


‘^cyl(j) 


‘^cyl(k) 







(E6) 
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For a cylinder the modified adjusted unit count such as <;^’cyl(i) distributed 

to the various columns for 1 = 0 to I = n, where n is the highest integral value less 
than the maximum normal component of impact velocity possible. Such procedure was 
necessary because oblique impacts would allow any normal velocity component from zero 
to the full value of s-s the case might be. That the values )^|> 

I v^l , and |vj^ | are the maximum possible normal components may be easily seen. 

For example, the surface of an i oriented cylinder would encoimter normal impact ve- 
locities dependent only on the vector velocity v^. The only component of v^ missing 
from v^ is parallel to the axis of the cylinder and hence can contribute nothing toward 
normal impact velocity on the surface. Hence, the maximum normal component of 
velocity on the surface would have the scalar magnitude | , the minimum would be 

zero. 

Sketch (d) illustrates the manner of partition of The i plane is shown, as 



i axis 
(d) 

well as its circular intersection with an i oriented cylinder of radius The angles 

at which the normal component of impact velocity on the cylindrical surface will be equal 
to I and I+l are assumed as 6j and 

Now the normal impact velocity with impact at the angular position 6j will be 

= 1^1 cos 


or, with equal to I, 
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6j = cos 


I 


(E8) 
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Of many hypothetical impacts of meteoroids, each with a given velocity component v^j^, 
the fraction that will have a component of velocity normal to the cylindrical surface 
within the range from I to I + 1 must be equal to the distance between the two vectors 
v^ as shown in the sketch divided by the distance Hence, such fraction is 


cPl = 


^cyl 


(E9) 


From equations (E8) and (E9), 


and the part of the modified adjusted unit count 
column of table VI(b) is 

•^cyKi)! = <Pcyl(i) <^1 

For any case in which I + 1 is greater than and I is smaller than | Vj^ | , zero must be 
substituted for the second term in the right-hand side of equation (ElO). When I is 
equal to or greater than | vT^| , <Pcyx(j)i must be arbitrarily zero without reference to 
equations (ElO) and (Ell). 

For the three sample meteors, the values of ‘>^cyl(j)’ ‘^cyl(k) 

with use of equations (E6) are shown in table V. Also shown, for each meteor, are the 
values of n = 0, 1, 24, 25, 26, and 79. Those values of ^cyl(i)n 

entered in the proper positions in the skeleton tabulation shown in table VI(b) . 

After such entries in the full tabulation were completed for all meteors treated, the 
procedure was exactly as with planes. For j and k orientations of cylinders, similar 
tabulations were made with use of equations analogous to (ElO) and (Ell). The results 
are tabulated in table IV and plotted as the circular points in the three parts of figure 2. 


1 - 


'I + 1 


(ElO) 


V, 


jk' 


‘^cyl(i) 


that must be entered in the 1 
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III 


Distribution of Normal Component of Impact Velocity Upon Spheres 

A tabulation not shown even in skeleton form but similar to those of table VI was 
provided for all impacts, normal or oblique, on the surface of a sphere. The results 
were as shown in the last column of table IV and as shown by the plotted circular points 
in figure 3. (The second column of table IV is a recapitulation of results explained and 
discussed in ref. 4.) 

As a sphere presents the same size and shape of surface with the same relative 
orientation for a meteoroid flux from any direction, no directional characteristics of the 
relative velocity Vq need be considered. Any meteor can impact a sphere with normal 
component of the impact velocity within the range from zero to the full value of v^. No 
consideration needs to be given to a modification of for variation in effective vul- 
nerable area similar to those that produced the modified adjusted unit counts 
^cyl(i)' *^w disable directly instead of or <Pcyi(j)- 

The value of the adjusted unit count (p^, however, does need to be distributed among 
the various values of I from 1 = 0 to I = n, where n is the largest integer that is less 
than Vq. Sketch (e) illustrates the principle of such distribution. In that sketch, a 
hemisphere is shown above an intersecting plane. Such plane may be arbitrarily selected, 
but containing an axis that passes through the center of the sphere in a direction parallel 
to Also shown, at the right, is a plane normal to the direction of Projection 

lines extending from the hemisphere to the plane at the right represent three concentric 
cylindrical surfaces: (1) circumscribing the hemisphere, that is, with all elements of 



(e) 
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the cylindrical surface tangent to the hemisphere, (2) with all cylindrical elements pass- 
ing through the circular locus on the hemisphere for values of the angle 6 equal to 
and (3) the same as (2) but for With angles 6j and having the same signi- 

ficance as earlier, equations (E7) and (E8) apply with substitution of Vq for |^^|- 
But, here, the probability that the normal impact velocity upon the sphere will be between 
I and I + 1 must be equal to the ratio of two areas within the plane at the right of the 
sketch. That ratio must be the area between the concentric cylinders that were enumer- 
ated (2) and (3), divided by the total area within the cylinder enumerated (1). Hence, 
equation (E9) is replaced with 


r2 


^sin^ 6j - sin^ ^I+l) 


(El 2) 


sph 


in which is radius of the sphere. Equation (ElO), then is replaced by 


<Pl 



I + 1\^ I I V 21+1 


'G 


(El 3) 


and equation (Ell) is replaced by the following expression for the fractional adjusted 
unit count for normal velocity components on a spherical surface, 


‘^sph(l) " '^w'^I 


(E14) 
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APPENDIX F 


CALCULATION OF CONCENTRATION FACTORS FOR METEOROID IMPACT 
UPON VARIOUSLY ORIENTED STRUCTURES 

For use in calculation of necessary armor thickness a concentration factor, or 
ratio, (pg is needed by which the expected frequency of impact of meteoroids upon one 
side of a randomly oriented surface may be multiplied to yield the frequency of impact 
upon a given structure with a given orientation. For a sphere, the ratio <p^ is always 
unity because in effect the sphere is a randomly oriented surface, having differential 
surface areas facing impartially in all possible directions. But, for planes and cylinders 
with i, j, and k orientations, the value of cp^ is never unity. 

For one example each of a plane and a cylinder the total of the next-to-last column 
in table V(a) or (b) provides an important portion of the desired ratio cp , which will be 
designated (pg. The significance of the factors signifies pi or cyl, and n 

signifies i, j, or k), as they appear in table V is explained in appendix E. The value of 
the factor <pg as defined is the total of those values of pertinent signifi- 

cance of m and n. From the explanation in appendix E, therefore, the value of (pg is 
the total of a modified adjusted unit count for all meteors treated, for the given orienta- 
tion of the structure, expressed as 


S ‘^m(n) ^ 

where the summation is performed for all meteors treated and where is a factor 

relating the probability of impact by the particular meteoroid under consideration upon 
the structure as oriented, relative to the probability of impact by that same meteoroid 
upon the same structure with an arbitrary basic orientation relative to the direction of 
motion of the meteoroid. For planes, that basic orientation was normal to the path of 
each meteoroid. For cylinders the basic orientation was with the cylinder axis parallel 
to any arbitrarily chosen line perpendicular to the meteoroid path. 

Now, from equation (FI), an average value of frequency of impact upon the structure 
as oriented, as a fraction or multiple of the frequency of impact with the basic orienta- 
tion relative to movement of each meteoroid, will be 

^or(av) “ ^ ^ 
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The manner of obtaining a value of from the value of <Por(av)’ planes and cyl- 
inders, will now be described. 


Concentration Factor for Planes 

For a plane ^ as used in equation (F2) will be regarded as the statistical fre- 
quency of impact per unit area on one arbitrarily chosen side of that area by all sporadic 
meteors from the single quadrant of space treated, in unspecified units, under the 
assumption that the chosen side of that area is always oriented to face each oncoming 
meteoroid normally. (With different units, ^(P-^ could as well be regarded as the fre- 
quency of impact upon one arbitrarily chosen side of a randomly oriented surface. ) 

Now the one side of a unit area arbitrarily oriented to face any oncoming meteoroid 
normally will receive the same frequency of impacts as the external surface of a sphere 
of one unit cross section (radius But the external surface of that sphere, in 

effect randomly oriented, will amount to four units of area. Hence, the frequency of 
impact upon it per unit area of outer surface will be one -quarter as great as for the 
chosen side of the always normally oriented unit area. So 

Z^w{r) = ^ E <^w 

where is the frequency of impacts upon one side of a randomly oriented surface 

by meteoroids from the single quadrant of space that was treated. By the same rea- 
soning. 


^w(r) ^ *^w 

where is the same adjusted unit count as <p^, but referred to one side of a ran- 

domly oriented area, instead of the always normal orientation. 

Now conversion of equation (F2) to an expression for frequency of impact relative 
to one side of a randomly oriented surface (for the one quadrant only) is analogous to a 
change of units. That is, equation (F2) will give the same result, but related to the 
random orientation rather than the normal orientation, if of equation (F4) is sub- 

stituted for cp^ in the denominator of the ri^t-hand side to ^ve 


'^g(qu) 


E‘^w(r) 


(F5) 
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Now of equation (F5) would be the desired end result (p^ except that only 

impacts from one quadrant of space have been considered. For one side only of the j 
or k oriented plane, adding the symmetrical flux from three additional quadrants of 
space would only double the frequency of impact, because meteoroids from two of those 
quadrants would impact the other side of each plane. But the impacts on a randomly 
oriented surface would be multiplied by four. Hence, for the j or k plane. 




*Pg(qu) 


‘P, 


g 


<P 


g 


2 E <^w(r) 1 £ 




w 


(F6) 


and, with the following definition. 


*i^pl(r) - - <^w 


(F7) 


equation (F6) becomes 


<Pg 


<Pl 


E ‘^pl(r) 


(F8) 


For either side of an i oriented plane, addition of symmetrical impacts from the 
four quadrants of space not treated would multiply the impacts by four, just as with the 
randomly oriented surface. Hence, for either side of the i plane, equation (F6) is re- 
placed by 




V. 


g 




g 


E ‘5^w(r) i ^ 


(p. 


w 


and equation (F7) becomes 


-^pKr) - 7 <Pw 


(F9) 


(FIO) 


for use in equation (F8). 

For the j plane as an example, values of from equation (F7), as shown in 
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Hence, the total of that column, 
Accordingly, as 

g‘ 


table V, were entered in the last column of table VI(a) 

3. 786x10'^, is the value of X/<Ppl(r) equation (F8) 

shown in the footnote of the table, equation (F8) yielded a value of 121. 07 percent for cp 
The procedure for the k plane was the same. For each side of the i plane the 
procedure was similar, but with use of equation (FIO) instead of (F7). Also, the sum- 
mation X/*^pl(r) *nade for all meteoroids treated, not just those that impacted the 
pertinent side of the i plane. 


Concentration Factor for Cylinders 

The procedure for an i, j, or k cylinder was like that for the j and k planes 
except that a different expression was required than equation (F7). For the i oriented 
cylinder, as an example, the value of cp' for use in equation (F2) was the total of the 
next-to-last column in table VT(b). 

With the basic orientation of a cylinder of length s-nd radius I'gyp the total 

cylindrical surface area would be 


^cyl ^’^^cyl^cyl 

and the equivalent area presented normally to a meteoroid path would be 

-^eff ^ ^’^cyl^cyl 

The same equivalent area would be offered by a sphere of radius 


_ -| /^^cyl^cyl 

"ph - y — 

The external surface of such sphere would have an area 

®sph “ '^’^^sph ~ ^^cyl^cyl 


(Fll) 


(F12) 


(F13) 


(F14) 


From equations (Fll) and (F14), 


= i (F15) 

Scyl 
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But the surfaces Sgp|j and should receive the same number of impacts, and the 

surface is in effect randomly oriented. Hence, the impacts per unit area upon the 

surface S^pj^ should exceed those upon S^yj by the ratio x/4. Accordingly, the 
equivalent of equation (F4) for use with the cylinders should be 


^■w(r) 4 (F16) 

Addition of impacts by meteoroids from three additional quadrants, for i, j, or k 
cylinders, would multiply total impacts by four, just as with a randomly oriented sur- 
face. Hence, by the same reasoning as with planes, equations (F7) and (F8) become 


"i^cyUr) 4 


(F17) 


and 




‘^cyl(r) 


(FI 8) 


The method of application of equations (F17) and (F18) in the example shown in 
tables V and VI(b) is identical to the example for a plane previously discussed. 
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o Adjusted and normalized count of meteors 



(d) Plane normal to apex of Earth movement, leading side. 
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(b) Plane normal to apex of Earth movement, trailing side. 

igure 1. - Distribution of norma! component of impact velocity of meteoroids upon a plane surface with 
three principal orientations. 



Adjusted and normalized count of meteors 
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(c) With line pointing toward Sun normal to plane. 


Adjusted and normalized count of meteors 



o 
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(d) Plane identical with ecliptic plane. 

Figure 1. - Concluded. 



o Adjusted and normalized count of meteors 

f lP"t''norm + 2.0)-2.0 

^ «v^orm) 0.0894 exp ^ 


8. 86x10-3 exp] 4 
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(a) With one end of cylinder axis pointing in direction of Earth movement. 


o Adjusted and normalized count of meteors 
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(b) Vi/ith one end of cylinder axis pointing toward Sun. 


o Adjusted and normalized count of meteors 

f 1 ["’"'norm * 2.0) - 2.050 f iF l" W' 
«*norm ' = 0. 0899 exp j 51555 J J + 7. 53x10 exp |-2 „ j 2 „ 
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(c) With cylinder axis normal to ecliptic. 

Distribution of normal component of Impact veixity of meteoroids upon cylindrical surface with three principal orientations 



o Adjusted and normalized count of meteors 
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Figure 3. - Distribution of normal component of impact velocity of meteoroids upon surface of a sphere. 
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Figure 5. - Product of concentration factor and average n^^ power of normal component of impact veixity of meteoroids on cylindrical 
surface. 


65 



(kmfsec) 



66 



NASA-Langley, 1970 31 E-5107 



rj cylinder 



Figure 8. - Effect of incidence angle of meteoroid impact on vulnerable area of plane or cylinder. 
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